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CHAPTER ONE 


Introduction 

Nature at the most fundamental level is described by four types of interactions, 
namely the gravitational, electromagnetic, weak and strong interactions. Gravity is 
best described by the theory of General Relativity, which was proposed in its hnal form 
a century ago on November 25, 1915 by Albert Einstein [1]. General Relativity is a 
classical held theory. The remaining three interactions are all described by quantum 
held theories. The theory of Quantum Electrodynamics was developed in the late 
1940s and early 1950s, largely by Sin-itiro Tomonaga, Julian Schwinger and Richard 
Feynman [2]. The theory of electro-weak interactions was put together in the 1960s 
by Sheldon Glashow [3], Steven Weinberg [4] and Abdus Salam [5]. The theory 
of strong interactions known as Quantum Ghromodynamics was discovered in 1973 
by David Gross and Frank Wilczek [6], and independently by David Politzer [7]. 
All these theories are experimentally well tested. But the unsatisfactory fact is that 
gravity is described by a classical theory while the rest of the fundamental interactions 
are described by quantum theories, hence the fundamental laws of nature do not 
constitute one logically consistent system. Despite the fact that this disparity does 
not affect physics as long as one stays safely below the Planck scale, one is forced 
into such a high energy regime when studying the very early universe or the final 
stages of star collapse. This requires a quantum held theory of gravity and perhaps a 
unihcation of all the fundamental interactions. The latter goal is pursued by String 
Theory, while there are quite a few different approaches to quantize gravity. 

1.1 Problems of Quantizing Gravity 

Soon after Quantum Field Theory was invented and applied to the electromagnetic 
held by Werner Heisenberg and Wolfgang Pauli [8], it was realized that perturbative 
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computation of the electromagnetic self-energy for an electron led to infinities. Pro¬ 
voked by Pauli, who thought the inhnities could be eliminated if gravitational effects 
were included, Leon Rosenfeld applied Quantum Field Theory to study the gravita¬ 
tional held in 1930 [9]. He showed that the gravitational self-energy of a photon is 
quadratically divergent in the lowest order of perturbation theory. Later in 1950, us¬ 
ing a manifestly gauge-invariant method due to Julian Schwinger, Bryce DeWitt [10] 
showed that this result implied merely a renormalization of electric charge rather than 
a non-zero mass of photon and it also implied that the one-loop contribution is zero. 

1.1.1 Unitarity 

In 1962, while studying the one-loop behavior of General Relativity quantized in 
a covariant gauge, Richard Feynman noticed that the diagrams are unitary only if 
hctitious particles are added besides the graviton [11]. They are hctitious because 
they were anticommuting bosonic particles circulating in the loop. Such nonphysical 
helds that are introduced in a manifestly Lorentz covariant quantization procedure 
to maintain unitarity at one-loop order of perturbation theory are usually referred 
to as “ghost” fields. DeWitt extended this to the two-loop case in 1964 [12] and to 
all orders in 1966 [13, 14, 15]. Ghosts were independently introduced by Ludwig D. 
Faddeev and Victor N. Popov in 1967 [16], who gave a prescription for quantizing 
Yang-Mills gauge theories in general gauges. 

1.1.2 Renormalizability 

Even early on, Heisenberg noticed that the gravitational coupling constant G has 
the dimension of a negative power of mass and consequently the divergences in the 
perturbation expansion of General Relativity will be different from renormalizable 
quantum held theories. 

In 1974, ’t Hooft and Veltman investigated all the one-loop divergences in Gen¬ 
eral Relativity. They discovered no physically relevant divergences in the one-loop 
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S-matrix. However, the finiteness of pure gravity was destroyed by coupling a single 
scalar to it [17]. Soon after that, essentially all of the known matter couplings were 
investigated and it was found that none of the matter fields would share the one-loop 
finiteness of pure gravity. In 1986, Marc Goroff and Augusto Sagnotti made it clear 
that our conventional quantization techniques do not work for gravity [18]. 

Non-renormalizability means simply that the theory is effective and its predictions 
can be trusted only at sufficiently low energies. For example, the Schrodinger equation 
is non-renormalizable, yet its predictions are experimentally found to be correct. The 
same is true for the Four-Fermi Theory of weak interactions. Also, theories that are 
perturbatively non-renormalizable can be non-perturbatively renormalizable. If there 
is any hope at all for the renormalizability of General Relativity, that should come 
from a non-perturbative approach. 

1.2 Semi-classical Approach 

At the one-loop level, the contribution of the gravitational loop is of the same 
order as the contributions of matter fields. At usual energies less than the 10^® GeV 
scale, the contributions of additional gravitational loops are highly suppressed. There¬ 
fore, a semi-classical concept applies when the quantum matter fields together with 
the linearized perturbations of the gravitational field interact with the background 
gravitational field [19, 20, 21]. 

Not withstanding the difficulties of quantizing gravity, Jacob Bekenstein and 
Stephen Hawking took up semi-classical studies of black holes in the early 1970’s. 
In such an approach, the gravitational field is treated as a classical background on 
which quantum fields are studied. It was shown that black holes behave like black 
bodies and emit thermal radiation [22, 23]. Hence one can assign temperature and 
entropy to a black hole. In fact, they showed entropy is proportional to the area of 
the event horizon [24, 25]. This started the modern era of quantum gravity. 
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The semi-classical theory of the early universe called Inflationary theory [26] has 
been very successful too. While black holes are the theoretical laboratories for quan¬ 
tum gravity, Inflation is the experimental laboratory that might give us the first 
glimpses of quantum gravity [27, 28]. The recent observational constraints [29] on the 
inflationary gravitational waves imprinted in the cosmic microwave background are 
strengthening our hopes that experimental quantum gravity is in sight. 

1.3 Common Ground 

In the absence of experimental data, it is pragmatic to see what common features 
are shared by different approaches to solve the problems in quantum gravity. Today, 
reproducing the black hole entropy formula has become a standard test for any theory 
of quantum gravity. String Theory has successfully reproduced these results, although 
only in certain extremal black hole cases [30]. In Loop Quantum Gravity, the black 
hole entropy has been computed for the Schwarzschild black hole, except for a factor 
of one fourth in the formula, which can be obtained with an appropriate choice of 
a free parameter of the theory [31]. It is often argued that another common feature 
shared by theories of quantum gravity is the similar fractal behavior [32]. It is shown 
that many, if not all, approaches to quantum gravity predict a spectral dimension 
of 2 in ultraviolet regime. However, it has recently been shown that Bosonic String 
Theory behaves differently [33]. One has to investigate such a spectral behavior of 
Superstring Theory before drawing any hnal conclusions. 

1.4 Current Status of Quantum Gravity 
To the best of our knowledge today, one cannot have a quantum held theoretic 
description of General Relativity without an idea from outside the conventional frame¬ 
work. One such idea that has been pursued for many decades is supersymmetry. In 
1983, it was established that A/” = 4 super-Yang-Mills theory in four dimensions is 
ultraviolet hnite to all-loop orders [34]. The supersymmetric version of General Rel- 


4 



ativity called supergravity showed improved ultraviolet behavior in the early days 
although it was thought to diverge in the third loop. Now, contrary to all the power¬ 
counting arguments, it is shown that J\f = 8 supergravity in four dimensions is hnite 
to seventh-loop order in perturbation theory [35]. That is a remarkable achievement. 

Another area of active research is String Theory whose most fundamental result 
to date is that under certain caveats, gravity is dual to certain conformal held theo¬ 
ries [36]. This result might contain some deeper clues to quantum gravity. 

The latest progress comes from a new way of describing gravity within the frame¬ 
work of quantum held theory by breaking Lorentz and diheomorphism symmetries. 
This was proposed by Petr Hofava in 2009 [37]. In the relatively short span of its ex¬ 
istence, it has been shown that this theory is a part of String Theory and at the same 
time makes contact with other seemingly diherent approaches to quantum gravity. 

Despite making enormous progress in our understanding of the high energy be¬ 
havior of gravity, our quest for quantum gravity is far from over. 
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CHAPTER TWO 


Gravitational Theories with Anisotropic Scaling 
2.1 Introduction 

In the previous chapter, we listed the core problems in quantizing gravity and 
argued that in order to achieve this feat, we need ideas from outside the conventional 
framework. The investigations in this dissertation are concerned with the latest of 
such ideas, inspired by anisotropic scaling in the theories of condensed matter physics. 
Before diving into the details of gravitational theories with anisotropic scaling, we 
discuss some important ingredients of the theory and why they are essential. 

2.1.1 Higher-order Terms 

The non-renormalizability of General Relativity means that it is an effective the¬ 
ory and the Einstein-Hilbert action contains only the terms relevant at low energies. 
One is naturally tempted to add higher-order curvature terms to the action, thereby 
making the theory applicable at high energies. This possibility was first explored in 
1962 by Ryoyu Utiyama and Bryce S. DeWitt [38]. They noticed that the action of 
quantum gravity should contain functionals of higher derivatives of the metric ten¬ 
sor besides the Einstein-Hilbert action. But is such a theory renormalizable? This 
question was answered affirmatively in 1977 by Kellogg Stelle [39]. He showed that 
the theory is renormalizable with quadratic curvature invariants. However due to the 
presence of higher time derivatives, such a theory has the negative norm state called 
“ghosts” which allow the probability to be negative and hence breaks unitarity. In 
fact, as far back as in 1850, Mikhail Ostrogradsky showed that the presence of time 
derivatives higher than two generally leads to the problem of ghosts [40]. Combining 
all these ideas, Hofava added only the terms containing higher spacial derivatives 
while keeping the time derivatives to the second order. Also, power-counting renor- 
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malizability restricts the number of spacial derivatives to at least six. This means that 
space and time are treated differently and hence the local Lorentz and diffeomorphism 
symmetries are broken. 

2.1.2 Lorentz Invariance 

Invariance under the Lorentz symmetry group is a cornerstone of modern physics, 
and is strongly supported by observations. In fact, all the experiments carried out so 
far are consistent with it [41], and no evidence shows that such a symmetry must be 
broken at certain energy scales, although it is arguable that such constraints in the 
gravitational sector are much weaker than those in the matter sector [42]. It should 
be emphasized that the breaking of Lorentz invariance can have signihcant effects on 
low-energy physics through the interactions between gravity and matter, no matter 
how high the scale of symmetry breaking is [43, 44]. The question of how to prevent 
the propagation of the Lorentz violations into the Standard Model of particle physics 
remains challenging [45]. Nevertheless, it is known that demanding renormalizability 
and unitarity of gravity often lead to the violation of Lorentz invariance [46]. There 
are benehts of violating Lorentz invariance [47]. Recently, a mechanism of super- 
symmetry breaking by coupling a Lorentz-invariant supersymmetric matter sector to 
non-supersymmetric gravitational interactions with Lifshitz scaling was shown to lead 
to a consistent HL gravity [48]. 

2.1.3 Diffeomorphism Invariance 

In many theories of quantum gravity and especially in String Theory, it is believed 
that the concept of space and time is an emergent notion that is not present in the 
fundamental theory [49]. Another way of looking at it is that one does not have to 
respect the symmetries of spacetime that are thought to be fundamental at known 
energies, such as the local Lorentz and diffeomorphism symmetries, especially given 
the fact that our understanding of spacetime at the Planck scale is highly limited. 
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2.1.4 Ultraviolet Fixed Point 

We know from quantum field theory that the value of coupling strengths of fun¬ 
damental interactions depend on the energy scale at which they are measured. This 
variation of coupling strengths with energy scale is called the running of the couplings. 
This has been experimentally seen for the Standard Model interactions. The essen¬ 
tial tool for studying this is renormalization group flow, developed by Kenneth G. 
Wilson [50], which is a sophisticated version of the general theory of thermodynamic 
phase transitions due to Lev Landau from 1937 [51]. This is closely related to criti¬ 
cal phenomena in statistical mechanics [52], originally developed by E. M. Lifshitz in 
1941 [53]. A quantum held theory is said to have an ultraviolet hxed point if its renor¬ 
malization group how approaches a hxed point at very high energies. Theories with 
the trivial or Gaussian ultraviolet hxed point are said to be asymptotic free theories. 
They may be, but not necessarily are, strongly-coupled theories whose study require 
nonperturbative methods. The famous example is QGD. Hofava-Lifshitz Theory is 
also thought to be of this kind. Theories with a nontrival or non-Gaussian hxed point 
are said to be asymptotic safe theories. A notable example is asymptotic safety in 
quantum gravity proposed by Steven Weinberg [54]. 

2.2 Hofava-Lifshitz Theory 

2.2.1 Original Version 

Hofava achieved all the above mentioned objectives in 2009 [37] by borrowing 
from condensed matter physics the idea of anisotropic scalings [53] of space and time, 
given by 

X —)■ f'x, t ^ Ft, (2.1) 

where 2 : is called the dynamical critical exponent. This was hrst used to construct 
scalar held theory by E. M. Lifshitz, hence the theory is often referred to as the 
Hofava-Lifshitz (HL) Theory. The anisotropic scaling provides a crucial mechanism; 



the gravitational action can be constructed in such a way that only higher-dimensional 
spatial, but not time, derivative operators are included, so that the UV behavior of 
the theory is dramatically improved. In (3 -|- l)-dimensional spacetimes, HL theory 
is power-counting renormalizable provided that z > 3 [37, 55, 56]. In this disserta¬ 
tion, we will assume that z = 3. At low energies, the theory is expected to flow to 
z = 1. In this limit the Lorentz invariance is “accidentally restored.” The exclusion 
of high-dimensional time derivative operators, on the other hand, prevents the ghost 
instability, whereby the unitarity of the theory is assured. 

The anisotropy between time and space mentioned above is conveniently expressed 
in terms of the Arnowitt-Deser-Misner (ADM) [57] variables: the lapse function N, 
the shift vector and the three-dimensional metric dehned on the leaves of constant 
time gij where i, j = 1, 2, 3. In the metric form, it is given by 

= -N'^c^df + gij {dx^ + N^dt) {dx^ + N^dt) , {i, j = 1,2,3). (2.2) 

Under the rescaling (2.1) with z = 3, N, and gij scale as 

N^N, g,,^g,,. (2.3) 

The requirement that the foliation dehned by these leaves be preserved by any 
gauge symmetry implies that the theory is covariant only under the action of the 
group Diff(M, J^) of foliation-preserving diffeomorphisms, 

6 t = -f{t), 5 x^ = -CM, (2.4) 

for which N, and gij transform as 

^9ij ^iCj T ^jCi T fdijt 

SN, = NkV.C’^ + C’^VkN. + gd’^ + NJ + Nj, 

6 N = + iV/+ iV/, (2.5) 

where / = df/dt, Vj denotes the covariant derivative with respect to the 3-metric 

gij, and Ni = gi^N^, etc. From these expressions one can see that the lapse function 
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N and the shift vector iVj play the role of gauge fields of the Diff(M, J^) symmetry. 
Therefore, it is natural to assume that N and iVj inherit the same dependence on 
spacetime as the corresponding generators 

N = N{t), N, = N,{t,x), (2.6) 

which is clearly preserved by Diff(M, J^) and usually referred to as the projectability 
condition. 

The non-projectable version of HL gravity [58] is self-consistent and passes all the 
solar system, astrophysical and cosmological tests. In fact, in the infrared the theory 
can be identihed with the hyper surface-orthogonal Einstein-aether theory [59] in a 
particular gauge [60, 61, 62], whereby the consistence of the theory with observations 
can be deduced. 

2.2.1.1. Problems. Due to the restricted diffeomorphisms, one more degree of 
freedom appears in the gravitational sector, i.e., the spin-0 graviton. This is poten¬ 
tially dangerous, and needs to decouple in the infrared in order to be consistent with 
observations [63]. In particular, the spin-0 mode is not stable in the Minkowski back¬ 
ground in the original version of the HL theory [37] and in the Sotiriou, Visser and 
Weinfurtner (SVW) generalization [64, 65]. However, in the SVW setup, it is sta¬ 
ble in the de Sitter background [66]. In addition, non-perturbative analysis showed 
that it indeed decouples in the vacuum spherical static [67] and cosmological [68] 
spacetimes. 

Another potential complication of the HL theory is that the theory becomes 
strongly coupled when the energy is very low [69]. However, as long as the the¬ 
ory is consistent with observations when the nonlinear effects are taken into account, 
this is not necessarily a problem, at least not classically. A careful analysis shows 
that the theory is consistent with observations in the vacuum spherically symmetry 

static case [67] and in the cosmological setting [70, 71, 72]. 
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2.2.2 Extended Symmetry 

One way to overcome the above problems is to introduce an extra local U{1) 
symmetry, so that the total symmetry of the theory is extended to [73] 

U{1) K Diff(M, T). (2.7) 

This is achieved by introducing a gauge held A and a Newtonian prepotential cp. 
One consequence of the U{1) symmetry is that the spin-0 gravitons are eliminated 
[73, 68]. Effectively, the spatial diffeomorphism symmetries of General Relativity are 
kept intact, but the time reparametrization symmetry is contracted to a local gauge 
symmetry [74]. The restoration of general covariance, characterized by Eq. (2.7), 
nicely maintains the special status of time, so that the anisotropic scaling (2.1) with 
z > 1 can still be realized. As a result, all problems related to them, such as instability, 
strong coupling, and different propagation speeds in the gravitational sector, are 
resolved. 

Under the Diff(M, T), the helds A and p transform as 

M = C%A + fA + fA, 

Sp = fp + Cdip, (2.8) 


while under the local U(l), the gauge helds together with N, Ni and gij transform as 


6aA 

= a-N^Via 

6aP 

= —a, 

5aN 

= 0, 


= NVia, 

daQij 

= 0, 


(2.9) 


where a is the generator of the local U{1) gauge symmetry. 


A remarkable by-product of this “non-relativistic general covariant” setup is that 

it forces the coupling constant A, introduced originally to characterize the deviation 
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of the kinetic part of the action from General Relativity [37], to take exactly its 
relativistic value A = 1. A different view can be found in [75]. The 77(1) symmetry 
was initially introduced in the case of A = 1, but the formalism was soon extended 
to the case of any A [75, 76, 77]. In the presence of a 77(1) symmetry, the consistency 
of the HL theory with solar system tests and cosmology was systematically studied 
in [78, 79, 80]. In particular, it was shown in [81] that, in order for the theory to be 
consistent with solar system tests, the gauge field A and the Newtonian prepotential 
(f must be part of the metric in the infrared limit. This ensures that the line element 
ds^ is a scalar not only under Diff(M, but also under the local 77(1) symmetry. 
In this model, the formation of black holes by gravitational collapse was hrst studied 
in [82]. 

Another possibility is to give up the projectability condition (2.6). This allows 
for new operators to be included in the action, in particular, operators involving 
Oj = N^i/N. In this way, all the problems mentioned above can be avoided by properly 
choosing the coupling constants. However, since this leads to a theory with more than 
70 independent coupling constants [83], it makes the predictive power of the theory 
questionable, although only five coupling constants are relevant in the infrared. 

A non-trivial generalization of the enlarged symmetry (2.7) to the nonprojectable, 
N = N(t,x), case was discovered by Zhu, Wu, Wang and Shu [84, 85]. It was shown 
that, as in General Relativity, the only degree of freedom of the model in the gravita¬ 
tional sector is the spin-2 massless graviton. Moreover, thanks to the elimination of 
the spin-0 gravitons, the physically viable range for the coupling constants is consid¬ 
erably enlarged, in comparison with the original non-projectable version, where the 
extra 77(1) symmetry is absent. Furthermore, the number of independent coupling 
constants is dramatically reduced from more than 70 to 15. The consistency of the 
model with cosmology was recently established in [85, 86, 87]. In the case with spher¬ 
ical symmetry, the model was shown to be consistent with solar system tests [88]. 
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In contrast to the projectable case, the consistency can be achieved without taking 
the gauge held A and Newtonian prepotential to be part of the metric. Finally, the 
duality between this version of the HL theory and a non-relativistic quantum held 
theory was analyzed in [89], and its embedding in String Theory was constructed 
in [90]. For other examples see [91]. 

It is remarkable to note that, despite of the stringent observational constraints 
on the Lorentz violation [41], the non-relativistic general covariant Hofava-Lifshitz 
gravity constructed in [84] is consistent with all the solar system tests [92, 93] and 
cosmology [80, 94]. 


Table 2.1. A summary of various versions of the Hofava-Lifshitz Theory 


Property of the theory 

Diff(M, H) 

U{1) X 

Diff(M, H) 

N = N{t) 

N = N{t,x) 

N = N{t) 

N = N{hx) 

Independent Couplings 

11 

>75 

11 

15 

No Spin-0 Graviton 

X 

X 

/ 

/ 

No instability of M 4 

X 

/ 

/ 

/ 

No Strong Coupling 

X 

/ 

/ 

/ 

Solar System 

/ 

/ 

/ 

/ 

Cosmology 

/ 

/ 

7 

/ 

String Theory Embedding 

? 

? 

7 

/ 


2.3 Non-relativistic General Covariant HL Theory 
In this section, we shall give a very brief introduction to the non-relativistic general 
covariant theory of gravity with projectability for an arbitrary coupling constant A 
and the enlarged symmetry (2.7). All our discussions in this dissertation are based 
on this version of the theory. 

The total action is given by 

S = J dt X N — C,y -\- C,^p -|- L C,\ C,m) •, (2.10) 
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where g = detgfjj, and 


Ck = 


Ca 

Cx 


ipQ''^ ( 2Kij + ViVjip ), 

A 


N 


2A,-R), 


(l-A) (vV)V2irvV 


( 2 . 11 ) 


Here the coupling constant Ag, which acts like a three-dimensional cosmological con¬ 
stant, has the dimension of (length)”^. The Ricci and Riemann terms all refer to 
the three-metric Qij. Kij is the extrinsic curvature, and Qij is the 3-dimensional 
“generalized” Einstein tensor dehned by 

Rij ~ 7^ i~9ij + VjlVj -|- VjiVj), 

1 

Qij = Rij ——QijR + AgQij. ( 2 - 12 ) 


Cm is the matter Lagrangian density and Cy denotes the potential part of the action 
given by 

Cv = C^fl'o + giR + ^ (fi'2-R^ + gsRijR^^) 

+1 {g,R^ + g,R + geR^RiR’l) 

[gjRV^R + gs , (2.13) 

which preserves the parity, where the coupling constants (s = 0,1, 2 ,..., 8) are all 
dimensionless. The relativistic limit in the IR requires that 


gi = -C C 


1 

levrG’ 


( 2 . 14 ) 


where G denotes the Newtonian constant. 

Variation of the total action (2.10) with respect to the lapse function N(t) yields 
the Hamiltonian constraint 


Ck + Cv- (pQ^^V^Vgip - (1 - A) (VV)' 
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SttG 


d^x^/g.R, (2-15) 



where 


rt JiNCM) 


r = 2 - 


SN 


(2.16) 


Variation of the action with respect to the shift iV® yields the super-momentum 
constraint 


7^ij - ‘fGij - (l - = STiGJi, 


where the super-momentum vr®-^ and matter current J® are dehned as 

7r®t = -K^^ + XKg^\ 

J‘ ^ -/vhli. 

SNi 

Similarly, variations of the action with respect to cp and A yield, respectively. 


( 2 . 17 ) 


( 2 . 18 ) 


e«(A'i,+ V,V,h + (l-A)vhA' + Vy) = SwGJ^, (2.19) 

R-2A, = SnGJA, (2.20) 


where 


j _ SCm j _ JiNCM) 

^ to - 1 ^ A - " " 


( 2 . 21 ) 


Sip ’ SA 

On the other hand, variation with respect to Qij leads to the dynamical equations 

- (1 - I = -2 {kY + 2\KK^^ 

+ ^Vfc [N^J - -2{l- \)\{K + VV) V*VV + 


+ (1-A) 

+ + fH + SttGJ, 

where = {J + fj /2, and 


+ 2 (+ C ,\) g''^ 


( 2 . 22 ) 


p®j _ 


1 - 
iV 


- (V*V^' - g^^y^]A 


pij = P*-?' 


n=l 




Jg 5g. 


(2.23) 




s=0 
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with Us = (2, 0, —2, —2, —4, —4, —4, —4, —4). The 3-tensors and are given 
below. 

(-Fo)ij = 

/ , 1 

(-Fi)ij = ——QijR -|- Rij, 

{^2)ij = + ‘^RRij — "2^/ {jS j)R + ‘^dij'^'^R, 

{Fsh = -^g^JRmnR^" + 2R,kR■-2V’^Vi^RJ)k + V^R^J+g^JVmVnR^^, 

{Fi)ij = ——gijR^ + 3R^Rij — 3'’^(^iVj)R^ + 3gij'’^‘^R^, 

(Fs)., = + + 

-2V"V(*FF,> + (FF,,) + (FF™”), 

(Fe),, = ~g,,R^R^^Rl, + 3F™"F,,F^, + (F,„F”) + (F^F'") 

-3VfcV(, (F,>F"'=) , 

(F?)*^ = —gij(VRY — (VjF) (VjF) -|- 2FjjV^F — 2V(jVj)V^F -I- 2gij'V‘^R, 

{Fs)^, = -^g^, {VpR^n) - V"F,, + (V^F^O (V,F”^") 

+2 (VpF,„) (V^F;^) + 2V"V(,V2F,)„ + 2V„ {RlV^^Rf^) 

-2Vn (F^(,V,)F”^") - 2Vn (F^(,V”F”;) - g,,V^V^V^Rmn, (2.24) 





where 


F I (2A- + V^tj) B« - 2 (2A'f + VWiip) fl“ - 2 ( 2 AI + V'Vtfj) W'‘ 

- (2A5 - fl) ('2A-‘^ + V‘VV) } , 




{2ViV'7j'‘ - v7? - 


+ v6"'(^ + vv) 
(VWi/i'Os'-'}. 


-■pe”(—+ vV 




(2.25) 


/^' = ^|(2F*^' + V*VV) -^(2F + VV)l?*^j. (2.26) 
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The stress 3-tensor is defined as 


= 


2 5{^Cm) 

y/9 ^9ij 


(2.27) 


The matter quantities (J*, J\ Ja, t*-^) satisfy the conservation laws 


d^Xy/g 




= 0 , 

(2.28) 




1 Tfc yy j 

-fn^ (-/gJi),. - - v.ivj - jiv.X + - ^v,A = 0, 

(2.29) 


as 


In General Relativity, the four-dimensional energy-momentum tensor is defined 

1 5 




dg, 


(4) 

fiiy 


where /i, z/ = 0,1, 2, 3, and 


Introducing the normal vector to the hypersurface t = constant by 




one can decompose as follows [95]: 


Ph = Si = 

g.. = T 

^IJ - fll /' ^2 ’ 


(2.30) 


gi"J = -N^ + N^N,, gl,f = N,, g^=g,y (2.31) 


(2.32) 


(2.33) 


where is the projection operator defined by = g^^} +n^nu. In the relativistic 
limit, one may make the following identification: 


("-f ) Jii ^ij'} ( ‘^PHi ^ij) ■ 


(2.34) 
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2-4 Summary 

In this chapter, we discussed how Hofava-Lifshitz Theory is constructed and the 
reasons for including certain features into the theory. After outlining the aspects of 
various versions of the theory, we have given a brief introduction to the projectable 
version with extended symmetry for an arbitrary value of the coupling constant. The 
rest of the dissertation deals only with this version of the theory. 
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CHAPTER THREE 


Spherically Symmetric Spacetimes 

This chapter published as: J. Greenwald, V. H. Satheeshkumar and A. Wang, 
“Black holes, compact objects and solar system tests in non-relativistic general 
covariant theory of gravity,” JCAP 1012, 007 (2010). 

3.1 Introduction 

In this chapter, we investigate systematically the spherically symmetric spacetimes 
and develop the general formulas of such spacetimes hlled with an anisotropic fluid 
with heat flow. We study the vacuum solutions. Although the solutions are not 
unique, when we apply them to the solar system, we find that these tests pick up the 
Schwarzschild solution generically. Later, we consider the junction conditions across 
the surface of a compact object with the minimal requirement that the matching is 
mathematically meaningful. In this chapter we work with A = 1. 

3.2 Spherically Symmetric Static Spacetimes 
The metric for spherically symmetric static spacetimes that obey anisotropic scal¬ 
ing and the projectability condition can be cast in the form 

ds'^ = —c^dt^ -I- [dr + e^~''dtY + r'^d'^Q, (3.1) 

in the spherical coordinates x* = (r, 9, 0), where = dO"^ + sin^ 9d(()^, and 

p = p{r), z/ = z/(r), W = {e^"^0,0}. (3.2) 
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The corresponding timelike Killing vector is ^ = dt. For the above metric, one hnds 





2z/' , , 

Rij 

= —^1^7+^- 

I^K 


Cu. 



rv' — (l — 


Vlij , 


Ca = 


2A 

3 


e^‘' (AgV^ - 1) + 1 (^(f" - z/V' + 2e^+V 

- 2 (^z/' - + 26^+^^) |, 

(1 - 2rz/') + (A^r^ - l) 




(3.3) 


s=0 


where a prime denotes the ordinary derivative with respect to its indicated argument, 
Qij = 6fSj + sin^ OSfSj, and £y^’s are given by 


£'”> = 2A 


2 e 


— 2u r 


rW _ 

— 


^( 2 ) _ 
— 


2e 


-Au 


‘^92 

4:94 


r(A _ 

— 


l^2^4 

2e-^^ 

^4^6 

2 e-®'^ 


2rz/' - (1 - 
2rv' - (l - 

2rv' - (l - 


+ 93 


+ 2(75 


3rV'2 - 2r(l - e^'^)z/' + (l - 

Qr^iy'3 _ 7^2^^ _ g2!.^^/2 + 4^('i _ 


+ 96 


5rV'^ - 3r2(l - + 3r(l - - (l - e^'') 


„2u\3 


(^4^6 


497 


2rV'(z/"' - 7z/'z/" + 6z/'3) - (l - e^’^)u"' - (9 - 7e^’^)u'u- 


+ 2{5-3e^'') + 


/3 


(1 - e^^)u" + 4z/' 2 + r(l - e^") V' + (l - 


-,2u\‘^ 


+93 







3r^ 

’(z/"-4z/'2)z/" + 4z/'^' 

-2r3(z/"-2z/'2)z/' + r2 

4(l-e2-)z/"- 

(3-8e2^)z/'2 







+ 3r{l-e^'')i2' + Q{l-e^''y 


(3.4) 
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The Hamiltonian constraint (2.15) reads 


K + Cy — — SttGJ*) e''r‘^dr = 0, 


where 


{Agr"^ - l) + 1 {(p" - u'lp') -2{u' - AgVe^'') p! >, (3.6) 


while the momentum constraint (2.17) yields 


2rv' + (A^r^ - l) + 1 <p' = -87rGr^e^^’'-f^'^v, (3.7) 

where J* = (u, 0, O). It can be also shown that Eqs. (2.19) and (2.20) now 


{Agr^ - 1) + 1 (^(p" - z/V' + e^+>') - 2 (^z/' - 

= STrGr^e^V^, (3.8) 


2rz/' - {AgP - l) + ij = dTrGr^e^V^. 


The dynamical equations (2.22), on the other hand, yield 

2(^p' + z/^) + - + — {Cg, + Ca) 

= -re-2^ (Prr + + F4 + STTGe^'i 


(3.10) 


where 


p” + (2p' — z/') A ^a) 




Fee + Fgg + F^ + SirGr'^pe , 


(3.11) 


Tij = e^’^prSlSj + r'^pe^ij, 


F^ = 
o 


(A' + AF)Sl6^j + ^ p(A" - z/'A') + r(A' + Az/') - A(^1 - , (3.12) 





2 ^ 2 ^ 


and the -Fj/s and are given below. 

(fi).; = 

= - re-^''Fa,„ 

4r"(2i/"-3i.'") + (l-e2"')(7 + e2"') 

4r3(2z/"' - 7z/V" + 6z/'^) - 2rz/'(7 - - (l - (7 + e^'^) 


^—2v 


515] 


2e 


-4!3 


r 


(^ 3 ),, = ^[3r^(2z."-3z.'2) + (l-e^^)(5 + e^^) 


rr 
0 ^ Oj 


(^ 4 ),, = 


^ J.2 

4e-4^ 


r 


6 


3r^(z/"' - 7z/'z/" + 6z/'^) - 2rz/'(5 - 26^'^) - (l - (5 + 

16r^z/'(3z/" - 5z/'^) - 12r(l - (2rz/" - - 4z/') 




- (1 - e^'') (23 - 2262^^ - 


515] 


4e 


— 6 ;^ 


24r" 


_ 111,'^y + 10z/'4 


- 4r3(l7 - lSe^'')v'^ - 12r2(l5 - 

- (1 - e^'^) 12r^ [v'” - 7z/'z/") - 48rV" + 3r (^1 + 76^^^) v' 


-2(l-e2'^)(23 + e2^) 


rijj, 


(^ 5 ),, = 


2e 


-Au 


12 r 


4r'' 


v'y - llv'v” + 10z/'3) + v" 

(1 _ e^-) {y'" - 7z/V" + 6 z/' 3) - z/'(3z/" - 2//'^) 

+ r2(l - (2z/" - 15 z/' 2) + 4r(l - e^'') (l3 - 2e^y 

+ (l-e2")'(23-e2")L[5; 


2e 


— 61 ^ 


ISr" Iz/V'" + z/"2 - z/'2(llz/" - 10z/'2) 

^ 1 7(1 - z/'" - (53 - 496^^) z/V" + (45 - 42e^'^) v'^ 

24(1 - e^y - (97 - me^'')u‘ 

+ r(l - e^'') (13 - 15e‘^’')u' + 2(l - e2'^)^(l3 + e^’') }>%, 


/2 
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(Fe).. = ^|l0r3(3z/"-5z/'2)z/'-3r2 2(l-e2'^)z/" 

+ 3e2"z/'2] + 12r(l - - (l - e^'^)\u + e^") 

+ ^ 115^4 ^ 

- 3(1 - e^'')v''' - 3(1 - 7e^'')v'v'' - (25 + ISe^'^)//'^ 
+ 3r2 4(1 - e^'')v" - (l2 - lle^'^)//'^ 

+ 12r(l - e^'') (2 - e^'')F + 2(l - e2'^)^(l4 + |%, 

Sp-4!^ f r 1 

(F^)^. = J ^4 ^ _ 2^(4) ^ 20z/V"' + (I5z/" - 82 z/'2)z/" + 40//'^ 

+ 2r2 2(3 - - 3(5 - 

- 8r(3 - e"")//' - (l - e"") (7 + e^") |5[5; 

Qp-Gl' f r 

+ —^< + 16z/V(^) + (25z/" - 101z/'2)z/'" 

- (I27z/" - 326z/'2)z/V" - 120z/'® 

+ 2r^ (3 - e2^)z/"' - (33 - 7e^‘')u'u" + (45 - 6e^‘')u'^ 

- 2r^ 4(3 - e^'')u” - (51 - lle^^)z/'^ 

+ r(57 - 246^^ - e^'')F + 2(l - e^’') (7 + 


23 



(Fg) ^ - 68z/'z/"' - (59z/" - 358z/'=^) z/" - 224z/'^ 

+ 2r3(l3z/" - 29 z/'2)z/' - 8(5 - 2e^''y - 7(l3 - Ae^^)v'‘^ 

+ 16r(4 - e^^)v' + 6(l - (l + 

f r 

+ -^\ - (25z/" - lOlz/'^)//'" 

+ (I27z/" - 326 z/'2)z/'z/" + 120z/'^ 

+ r" [z/V'" - (5z/" - 13 z/'2)z/" - 14z/'"' 

-r^ 2(7 - e2'^)z/"' - 2(78 - 7e2'^)z/'z/" + 2(l07 - 6e2'^)z/'3 
+ 8(7 - e^'')u'' - (277 - 

- 16r(l3 - - 6(l - e^'^) (ll - |%, (3.13) 

= ^{[.^^(l-V^)-(l + r.')>" 

+ - 1 j + (3 + rv') v'ip' 

- 2e^+" [e^\^gr^ - l) + + 2e^+" (2 - r/i') z/'|5[(5; 

+ ^996“^'^I rv' — ip" + ^3 — rv' — e^"'^v'ip' — 2Agre‘^"(p' 

- 2e^+" (2 - rv' - + 4e^+" (v' - 1%, 
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ij = ^ I (^Agr^ - l) + 1 + (^ip' + (^z/'j p>' + 26^^+^ [ip' + ipp!^ 

+ Aipe'^+''[u' - A^re^^) - 2Agrpe^'^p>^5l5] 

+ ^e“^'^<|r99^z/' — Agre^'^^ip" + rp>[p' + + r{v' — Agre^'^^ip''^ 

- np {^ip' + z/'2 -(ip- 2re'^+^ - A^r z^V' 

— 2e'^^'' (a gr‘^e^''p>' + ^^z/'j + 2ripe'^~'^'^ (u' — Agre^'^j/i'lf2jj, 


[^iv, 3 )] ij = npu'ip" + ip'‘^ - ip (rz/' + 2) u'ip' + 2e^+'' (ip' + npv'p - 2pv''^ 515^ 

+ ]^e-^''\^r(p' - pv' + - rp(p' + 2e^+'^)z/" + 4r((^z/' - (^') + <^ z/V' 

_l_ gM+i" /irpl^i/' — p')i/' + 2pi/' + 2r{^p' — 3i/')p' 

Here we define a fluid with = pe as a perfect fluid, which in general conducts heat 
flow along the radial direction [96]. 

Since the spacetime is static, one can see that now the energy conservation law 
(2.28) is satished identically, while the momentum conservation law (2.29) yields 

vp' - {v' -p'r) -‘^{v-Pr+ Pe) + - \ Ja^ = 0- (3.14) 

To relate the quantities J*, J* and r^- to the ones often used in General Relativity, 
one can hrst introduce the unit normal vector to the hypersurfaces t = constant, 
and then the spacelike unit vectors and 0^, are defined as [95] 

rig = 5% 

x, = e% + e''5l, 

0^ = rsin0(5^. (3.15) 
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In terms of these four unit vectors, the energy-momentum tensor for an anisotropic 
fluid with heat flow can be written as 

= PHn^,ny +q[n^Xu + nuX^l)+PrX^lXu+Pe{0^9y + (t)^(t)y) (3.16) 

where pn, q, Pr and pe denote, respectively, the energy density, the heat flow along 
the radial direction and the radial and tangential pressures, measured by the observer 
with the four-velocity n^. Then, one can see that such a decomposition is consistent 
with the quantities J* and J\ dehned by 

Ph = —2J*, V = e^q. (3-17) 

It should be noted that the dehnitions of the energy density ph, the radial pressure 
Pr and the heat flow q are different from the ones dehned in a comoving frame in 
General Relativity. 

Finally, we note that in writing the above equations, we leave the choice of the 
17(1) gauge open. From Eq. (2.9) one can see that it can be used to set one (and only 
one) of the three functions A, cp and i\7 to zero. To compare our results with the one 
obtained in [73], without loss of the generality, we shall choose the gauge 

ip = 0. (3.18) 

Then, we hnd that 

^^ = 0. = = (3-19) 

3.3 Vacuum Solutions 

In the vacuum case, we have .P = v = Pr = pe = Ja = Jip = ^- With the gauge 
(3.18), from the momentum constraint (3.7) we immediately obtain v = constant, 
while Eq. (3.9) further requires 

zy = 0, Ag = 0. (3.20) 
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This is different from the solutions presented in [73], where u ^ 0, fi = —oo. Inserting 
the above into Eq. (3.8), it can be shown that it is satished identically. Since z/ = 0, 
from the expressions of {Fs)^p we find that = 0 for s 7 ^ 0 , and {Fq)^ - = —gij/2, 

so that 

Fij = -Agij, (3.21) 


where A = (‘^goj2. Substituting Eqs. (3.12) and (3.19)-(3.21) into Eqs. (3.8) and 
(3.9), we hnd that 


(2rfj,' + ije^^ = Ar^ — 2rA', 


/i" + 2/i'U' + - = 


3-2ai 


Ar — {rAy 


(3.22) 

(3.23) 


It can be shown that Eq. (3.23) is not independent, and can be obtained from 
Eq. (3.22). Therefore, the solutions are not uniquely determined, since now we have 
only one equation, (3.22), for two unknowns and A. In particular, for any given A, 
from Eq. (3.22) we find that 




2m 1,2 X 2 r . 

-1—Ar^ — 2A(r) H— / A(r )dr 

r 3 r 


(3.24) 


On the other hand, we also have 


4A' 

Ck = - 2 A, Cv = 2 A. (3.25) 

r 

Inserting it into the Hamiltonian constraint (3.5), we hnd that 

rA'{r)dr = 0. (3.26) 

Therefore, for any given function A, subject to the above constraint, the solutions 
given by Eqs. (3.20) and (3.24) represent the vacuum solutions of the HL theory. 
Thus, in contrast to General Relativity, the vacuum solutions in the HMT setup are 
not unique. 

When A is a constant (without loss of generality, we can set A = 0), from the 
above we hnd that 

1, / 2m 1 , . X 

A = 2 (V ^ 3 '^'’ ) ’ = 
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(3.27) 



which is exactly the Schwarzschild (anti-) de Sitter solution, written in the Gullstrand- 
Painleve coordinates [97]. It is interesting to note that when m = 0 we must assume 
that A > 0, in order to have /i real. That is, the anti-de Sitter solution cannot be 
written in the static ADM form (5.5). 


3-4 Solar System Tests 

The solar system tests are usually written in terms of the Eddington parameters, 
by following the so-called “parametrized post-Newtonian” (PPN) approach, intro¬ 
duced initially by Eddington [98]. The gravitational field, produced by a pointlike 
and motionless particle with mass M, is often described by the form of metric [99], 

= —e^'^c^dt^ -f e^^dr"^ + r‘^d^Vt^ (3.28) 


where T and $ are functions of the dimensionless quantity y = GM/ (rc^) only. For 
the solar system, we have GMq/c^ ~ 1.5 km, so that in most cases we have y 1. 
Expanding T and $ in terms of y, we have [99] 


= 



-I- ..., 


= 1 + 2 ^ ( 1 ^) +( 3 - 29 ) 

where (3 and 7 are the Eddington parameters. General Relativity predicts /5 = 1 = 7 
strictly, while the current radar ranging of the Gassini probe [100], and the procession 
of lunar laser ranging data [ 101 ] yield, respectively, the bounds [ 102 ] 


7-1 = (2.1±2.3) X 10■^ 

13-1 = (1.2± 1.1) X 10■^ (3.30) 


which are consistent with the predictions of General Relativity. 

To apply the solar system tests to the HL theory, we need first to transfer the 
above bounds to the metric coefficients /i and v. The relations between ($, T) and 



(/i, u) have been worked out explicitly, and are given by 


or inversely, 


p=-ln c - e j , z/ = <|> + \k. 


(3.31) 


^In 


(3.32) 


Inserting Eq. (3.29) into Eq. (3.31), we find that 


1, 2 fGM\ , . fGMY 


, . fGM\ 

- (')'-i)(^j + 


(3.33) 


Comparing Eq. (3.33) with Eq. (3.24) for A = 0, we find that in order to be consistent 


with solar system tests, we must assume that 


A{r) = O 


(3.34) 


Together with the Hamiltonian constraint (3.26), we find that this is impossible unless 
A = 0. Therefore, although the vacuum solution in the HMT setup is not unique, 
the solar system tests seemingly require that it must be the Schwarzschild vacuum 
solution. 

It should be noted that by choosing A{r) in very particular forms, the condition 
A = 0 could be relaxed [79]. But, such chosen A is not analytic (in terms of the 
dimensionless quantity y), and it is not clear how to expand it in the form of (3.33). 
Thus, in this chapter we simply discard those possibilities. 


3.5 Perfect Fluid Solutions 


In this section, let us consider a perfect fluid without heat flow, that is. 


Pr=Pe=P, n = 0. 
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(3.35) 



Then, together with the gauge choice (3.18), from Eq. (3.7) we hnd that v = constant. 
However, to be matched with the vacuum solutions outside of the fluid, we must set 
this constant to zero. 


z/ = 0, (3.36) 

from which we immediately hnd that Rij = 0, Fij is still given by Eq. (3.13), and 

Ca = 2AgA, 

9 A' 

= —^W,+r(rA:)%,. (3.37) 

Inserting the above into Eqs. (3.8)-(3.12), we hnd that 

J ip 

Ja = 

(r/)^ + 2rH'+ — Ar^ = 

-r/" + f + {rA'y + KgvA — Ar = 

where / = From the last two equations, we hnd that 

r^f - 2/ = -2r^ . (3.42) 

On the other hand, the conservation law of momentum (3.14) now reduces to 

pA^A'^O, (3.43) 

which has the solution, 

where po is an integration constant. 

Substituting it into Eq. (3.40), and then taking a derivative of it, we hnd that 

the resulting equation is exactly given by Eq. (3.42). Thus, both Eqs. (3.42) and 

(3.41) are not independent, and can all be derived from Eqs. (3.40) and (3.44). Then, 

in the present case there are hve independent equations, the Hamiltonian constraint 
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A. 


47rG’ 

—SttGPp, 

—STiGrp, 




(3.38) 

(3.39) 

(3.40) 

(3.41) 



(4.4), and Eqs. (3.38), (3.39), (3.40) and (3.44). However, we have six unknowns, 
A, /i, p, J*, J^, Ja- Therefore, the problem now is not uniquely determined. As in 
the vacuum case, we can express all these quantities in terms of the gauge held A. In 
particular, substituting Eq. (3.44) into Eq. (3.40) and then integrating it, we obtain 

Then, we hnd that 

4A' 

C^ = -2 (A - SvrGpo), Cy = 2A. (3.46) 

r 

Inserting the above into Eq. (3.5), we hnd that it can be cast in the form, 

p(r)dr = 0, (3.47) 



2m 

r 


+ ^ ^A — STrCpo^r"^ — 2A+- j A{r')dr' 


(3.45) 


where 


r = 


1 


27rG' 


47rG'po + 


A'{r) p{r) 


(3.48) 


From the above one can see that once A is given, one can immediately obtain all 
the rest. By properly choosing it (and p(r)), it is not difficult to see that one can 
construct non-singular solutions representing stars made of a perfect huid without 
heat how. To see this explicitly, let us consider the following two particular cases. 


3.5.1 Ag = D 

When Ag = 0, we have 

Jg, = Ja = ^, P = Po, (3.49) 

while p and J* are still given by Eqs. (3.45) and (3.48), respectively. To have a 
physically acceptable model, we require that the huid be non-singular in the center. 
Since Rij = 0, one can see that any quantity built from the Riemann and Ricci tensor 
vanishes in the present case. Then, possible singularities can only come from the 
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kinetic part, Kij, where the very hrst quantity is 


K = = + 2) 

/‘im S /" \ 

= - ( - h (A — StiGpo) — rA' — 3A+- / Air')dr' ) . (3.50) 

r \ r ^ J / 

Assuming that near the center A is dominated by the term r", we find that K is 
non-singular only when 

m = 0, a > 2. (3.51) 

For such a function A, Eq. (3.48) shows that is non-singular, as long as p{r) ~ 

C>(r2). 

3.5.2 A = Ao 

When A is a constant, from Eq. (3.44) we can see that the pressure p is also a 
constant. Then, the integration of Eq. (3.40) yields, 

A = (3.52) 

Inserting it into Eq. (3.42) we find that it is satisfied identically, while the Hamiltonian 
constraint (3.5) can also be cast in the form of Eq. (3.47), but now with 

Thus, the solutions of Eqs. (3.36) and (3.52) represent a perfect fluid with a constant 

pressure, p = p{Aq), given by Eq. (3.44). In this case, it can be shown that K is free 

of any spacetime singularity at the center only when m = 0 . 

It should be noted that in [70] it was shown that non-singular static solutions 

of a perfect fluid without heat flow do not exist. Since their conclusions only come 

from the conservation law of momentum, one might expect that this is also true in 

the current setup. However, from Eq. (3.14) we can see that in the present case the 

conservation law contains two extra terms, J^p and Ja- Only when both of them 

vanish can one obtain the above conclusions. Since in general one can only choose 
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one of them to be zero by using the gauge freedom, it is expected that non-singular 
static stars can be constructed by properly choosing the gauge held A. 

It should be also noted that the arguments presented in [70] do not apply to 
the case where the pressure is a constant. Therefore, when p' = 0 non-singular stars 
without heat how can be also constructed in other versions of the HL theory, although 
when p' ^ 0 this is possible only in the HMT setup. In addition, the dehnitions of the 
quantities p, p and v (= qe^) used in this chapter are diherent from the ones usually 
used in General Relativity. 


3.6 Junction Conditions 

To consider the junction conditions across the hypersurface of a compact object, 
let us hrst divide the whole spacetime into three regions, and S, where V~ (R’’") 
denotes the internal (external) region of the star, and S is the surface of the star. 
Once the metric is cast in the form (3.1), the coordinates t and r are all uniquely 
defined, so that the coordinates dehned in V~^ and V~ must be the same, = 

= (f, r,9,(l)). Since the quadratic terms of the highest derivatives of the metric 
coefficients p and u are only terms of the forms, z/"^, z/'V'" and the minimal 
requirements for these two functions are that z/(r) and p{r) are respectively at least 
and across the surface S, and that they are at least and elsewhere. 
For detail, we refer readers to [103]. Similarly, the quadratic terms of the Newtonian 
pre-potential are only involved with the forms, and ipcp'. Therefore, the 

minimal requirement for p is to be at least across the surface S. On the other 
hand, the gauge held A and its derivatives all appear linearly. Thus, mathematically 
it can even be non-continuous across S. However, in this chapter we shall require 
that A be at least too across S. Elsewhere, A and are at least C^. Then, we 
can write A and cp in the form, 

E{r) = E^{r-)H (r — ro) + E~ {r) [1 — 7/ (r — tq)] , (3.54) 
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where E = (A, </?), tq is the radius of the star, and H{x) denotes the 
function, defined as 

{ 1 , a; > 0, 

0 , a; < 0. 

Since A and ip are continuous (C°) across r = ro, we must have 

hmit,.^^+E+(r) = limit^^,.-E-(r). 

Then, we hnd that 

E'{t) = E^(r), 

E"{r) = E^,(r) + [E']-<5(r-ro), 

where 


[E']- 

E^{r) 


limit — limits 


r^rA ,r 


E+H{r -To) + E- 1-H{ 


r - ro) 


Combining the above with Eq. (6.6), we find that 

Ck = Ca = C^, 

Cv = C^ + C^v"5{r-ro), 
= C^ + Cl^6{r-ro), 


where 


£ Im _ 

V = 


8 ^ 76 -®^ - 




(^4^3 
-iu r 


2 rz/'- (1 - [z/"]". 




+ (1 - b']. 


J= J^ +J^-5(r-ro), 
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Heaviside 


(3.55) 


(3.56) 


(3.57) 


(3.58) 


(3.59) 


(3.60) 


Setting 


(3.61) 



where J = {J*, v, J^}, and has support only on S, we hnd that the 
Hamiltonian constraint (3.5) can be written as 


[Ck + Cv- - SttGJ*) e'^r^dr 


= —(87rGJ*’^”^ + 4 ™-£(rj, 


(3.62) 


where 

/ D / nVQ — e noo \ 

I{r)dr = limite^o ( / I{r)dr + / I{r)dr j . 

\Jo J TQ+e J 

It should be noted that in writing Eq. (3.62), we had used the conversion 


(3.63) 


y/gd^xf{r)6{r - tq) = /(ro). 


(3.64) 


The momentum constraint (3.7) will take the same form in Regions V^, while on 


the surface S it yields 


= 0 . 


(3.65) 


That is, the surface does not support impulsive heat flow in the radial direction. 
Similarly, Eqs. (3.8) and (3.9) take the same forms in Regions V^. While on S they 
reduce, respectively, to 


Agr^e^’' + (1 - e^^)] [ip']~ = SiiGr^e^'', 


= 0, (r = ro). 


(3.66) 

(3.67) 


On the other hand, the dynamical equations (3.10) and (3.11) take the same forms 
in Regions and on the surface S they yield. 


+ (l - e^")] [ip'Y + - ^o) 


= -2rR F/-+ 87rGe>^ 


(3.68) 




g 2 (i/-At) 


F',’” + 8irGrVJ’”), 


(3.69) 



where are given below. 


P'r = 4g,(3!''-;(l-e^‘'))Kl--8gr([!'®l--10!^Vr) 

+ - 3V[!/"]“^ 6(x) - (Sg, - 3gs)|i'"]“i5'(i)|, 

Fi^ = (2^il±MrLllFrS(x) + ^{48g42..' - (l - e^)] [.'T 

+ 2^5 18rz/' - 7(1 - [F']~ + 3ge 5rz/' - (l - [u"]~ 

— 85^7 — 16r^z/'— r^(25 {z/"}"''— lOlz/'^) [z/"] — 2(3 — [z/"] 

+C /8 3r2 [z/(^)]" - 48rV' [z/(=^)]" - 3r2(25 {z/"}+ - lOlz/'^) [z/"]" 

+ (rz/'— 14 + 2e^'^) [z/"] 15 ( 0 ;) 


'7 - 3^8)^ 

^4g6i2 


16z/' [z/"' ■ 


( 8^7 - 3^8)r r 


From Eqs. (3.54) and (3.57) we find that given by Eq. (3.14) takes the 

form 

F(,.,»)],, = F(.,„)]„ + [hl’«)].F(’'P'™> 

where 


K™)],, = ^ ^^‘'(i-Av^)-(1+ '■>'') wrm 

rv' “ (l “ g^*^) [f'V 

+y4’g“‘“' (g'- + 26"+*') [y fiy, 
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K3)1« = 2(^'-v»''+e'‘+‘')i^r 

-<^(^(^' + 26^+^^) [z/']" Qij. (3.71) 

L = L^ + L^^, (3.72) 

with L = (pr, pe). We can see that in general takes the form, 

lim ^ ^ - To) + L(2)^"^5"(r - tq). (3.73) 

The above represents the general junction conditions of a spherical compact object 
made of a fluid with heat flow, in which a thin matter shell appears on S. 

In the following we shall consider the matching of the perfect fluid solutions. Since 

zy± = 0, ip^ = 0, (3.74) 

we immediately obtain 

R± = 0. £J = 2A±, 4” = 0, 

= -A±4, Flp = 0, 

C% = 4’" = 0. (4),^ = (4™).^. = 0. (3.75) 

Then, from Eqs. (3.65), (3.66), (3.67) and (3.68) we hnd 

^ Jl^m ^ Jijn ^ _ q 

That is, the radial pressure of the thin shell must vanish. This is similar to what 
happened in the relativistic case [96]. 

In the external region, V~^, the spacetime is vacuum, and the general solutions 
are given by Eq. (5.75) 

1 r 2m 1 2 r 

= -In -h-A+r^ — 2yl’''(r) + - / A^{r')dr' , (3.77) 

2 r 3 r J 
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for which we have 


4A~^ 

£J, = ^- 2 A^, £i = 0 . 


(3.78) 


In the internal region, V , two classes of solutions of perfect fluid without heat 
flow are found and given, respectively, by Eqs. (3.45) and (3.52), which can be written 


jj, = ^ In < ^ ^A_ — — 2A (r) + - / A {r')dr' 


for A = 0, and 


A = ^ In I ^ A- - Ag Ao - SttGp 


for A = Aq, where p = A Aq/ 2 + po- Then, we hnd that 


— 


2(87rGpo-A_)+A^, A; = 0, 
2(87rGp - A_ + A“Ao), A~ = Aq, 


= 2 K-A-{r) 


(3.79) 


(3.80) 


To further study the junction conditions, let us consider the two cases A^ = 0 and 
A~ = Aq separately. 


3.6.1 A- = 0 

In this case, the continuity conditions of p and A across S read, 

1 A A 3 r° A A. A , 47rG 3 

m + TrAArjj + / AA{r)dr = - 7 ;—Port 

b J 3 

A+(ro) = /l“(ro), 


(3.81) 


where AA = A_|_ — A_ and AA = A"*" — A . Then, the Hamiltonian constraint (3.62) 


becomes 


pro roo 1 

p(r)dr + / rA+(r)dr = 

) Jro ’ ‘2 


(3.82) 


while the dynamical equation (3.69) reduces to 


AA = —87iG(po + , 


(3.83) 



where = pf^^^5{r — vq) [cf. Eq. (3.73)]. 

When the matter thin shell does not exist, we must set = 0, and 

Eqs. (3.81)-(3.83) become the matching conditions for the constants A±, m, po and 
the functions A^{r) and p{r). 

3.6.2 7l = 7lo 

In this case, it can be shown that the continuity conditions for p and A become 

1 r° 1 

m+-AArl+ A+(r)dr = roAo --{A~Ao + 87rGp)rl, 

A+{ro) = Ao, (3.84) 

while the Hamiltonian constraint (3.62) reduces to 

pro poo 

/ p{r)dr + 4 / rA'j^{r)dr = 2GJ*’^"^. (3.85) 

Jo Jro 

The dynamical equation (3.69), on the other hand, yields 

AA = - A; Ao - 87rG (^p + . (3.86) 

In all the above cases, one can see that the matching is possible even without a 
thin matter shell on the surface of the star, = q = by properly choosing 

the free parameters. 


3 .7 Conclusions 

In this chapter, we have systematically studied spherically symmetric static space- 
times generally hlled with an anisotropic fluid with heat flow along the radial direc¬ 
tion. When the spacetimes are vacuum, we have found solutions, given explicitly 
by Eqs. (3.20) and (3.23), from which one can see that the solution is not unique 
because the gauge held A is undetermined. When A = 0, the solutions reduce to the 
Schwarzschild (anti-) de Sitter solution. We have also studied the solar system tests 
and found the constraint on the choice of A. 
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It should be noted that we have adopted a different point of view of the gauge 
held A in the IR limit than that adopted in [73]. In this chapter we have considered it 
as independent from the 4D metric although it interacts with them through the 
held equations. This is quite similar to the Brans-Dicke (BD) scalar held in the BD 
theory, where the scalar held represents a degree of freedom of gravity, is independent 
of the metric, and its ehects to the spacetime are only throngh the held eqnations 
[104]. On the contrary, in [73] the anthors considered the gauge held R as a part of 
the lapse function, gu — —{N — A)'^ in the IR limit. 

We have also investigated anisotropic hnids with heat how, and found perfect huid 
solntions, given by Eq. (3.45). By properly choosing the gange held R, the solutions 
can be free of spacetime singularities at the center. This is in contrast to other versions 
of the HL theory [70] due to the coupling of the huid with the gauge held. We then 
have considered two particular cases in which the pressnre is a constant, quite similar 
to the Schwarzschild perfect huid solution. In all these cases, the spacetimes are free 
of singularities at the center. 

For a compact object, the spacetime outside of it is vacuum and matching condi¬ 
tions are needed across the surface of the star. With the minimal requirement that 
the jnnctions be mathematically meaningful, we have worked out the general match¬ 
ing conditions, given by Eqs. (3.62) and (3.65)-(3.69), in which a thin matter shell in 
general appears on the surface of the star. Applying them to the perfect hnids, where 
the spacetime ontside is described by the vacuum solutions (3.23), we have found the 
matching conditions in terms of the free parameters of the solutions. When the thin 
shell is removed, these conditions can also be satished by properly choosing the free 
parameters. 

Finally, we note that da Silva argued, in the HMT setnp, that the conpling con¬ 
stant A can still be diherent from one [105]. If this is indeed the case, then one might 
be concerned with the strong coupling problem fonnd in other versions of the HL the- 
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ory [69, 68, 106]. However, since the spin-0 graviton is eliminated completely here, as 
shown explicitly in [73, 105, 107], this question is automatically solved in the HMT 
setup even with A 7^ 1. It should be noted that da Silva considered only perturbations 
of the case with the detailed balance condition, and found that the spin-0 mode is 
not propagating. It is not clear if it is also true for the case without detailed balance. 
The problem certainly deserves further investigations. 
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CHAPTER FOUR 


Gravitational Collapse 

This chapter published as: J. Greenwald, J. Lenells, V. H. Satheeshkumar and 
A. Wang, “Gravitational collapse in Hofava-Lifshitz Theory,” Phys. Rev. D 88, 

024044 (2013). 

The study of gravitational collapse provides useful insights into the final fate of 
a massive star [108]. Within the framework of General Relativity, the dynamical 
collapse of a homogeneous spherical dust cloud under its own gravity was first consid¬ 
ered by Datt [109] and Oppenheimer and Snyder [110]. It was shown that it always 
leads to the formation of singularities. However, in a theory of quantum gravity, it 
is expected that the formation of singularities in a gravitational collapse is prevented 
by short-distance quantum effects. 

In this chapter, we study gravitational collapse of a spherical star with a finite 
radius in the HL theory with the projectability condition, an arbitrary coupling con¬ 
stant A, and the extra U{1) symmetry. In General Relativity, there are two common 
approaches for such studies. One approach relies on Israel’s junction conditions [111], 
which are essentially obtained by using the Gauss and Codazzi equations. An ad¬ 
vantage of this method is that it can be applied to the case where the coordinate 
systems inside and outside a collapsing body are different. Although Israel’s method 
was initially developed only for non-null hypersurfaces, it was later generalized to the 
null hypersurface case [112]. For a recent review of this method, we refer to [113] 
and references therein. The other approach is originally due to Taub [114] and relies 
on distribution theory. In this approach, although the coordinate systems inside and 
outside the collapsing stars are taken to be the same, the null-hypersurface case can 
be easily included. Taub’s approach was widely used to study colliding gravitational 
waves and other related issues in General Relativity [115]. 
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We follow Taub’s approach, as it turns out to be more convenient when dealing 
with higher-order derivatives. Moreover, in contrast to the case of General Relativity, 
the foliation structure of the HL theory implies that the coordinate systems inside 
and outside of the collapsing star are unique. Thus, also from a technical point of 
view, Taub’s method seems a natural choice for the study of a collapsing star with a 
hnite radius in the HL theory. 

4-1 Spherical Spacetimes Filled with a fluid 
Spherically symmetric static spacetimes in the framework of the HL theory with 
U (1) symmetry with or without the projectability condition are studied systematically 
in [78, 79, 81, 88, 116, 117, 118]. In particular, the ADM variables for spherically 
symmetric spacetimes with the projectability condition take the forms 

N = 1, 

gijdx^dx^ = dr"^ + difl, (4.1) 

in the spherical coordinates x* = (r, 9, 0), where dfl"^ = d9^ -|-sin^ 9 dcffl. The diagonal 
case W = 0 corresponds to fi(t,r) = —oo. On the other hand, using the U{1) gauge 
freedom, without loss of generality, we set 

(^ = 0, (4.2) 
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which uniquely hxes the gauge. Then, we hud that 


= o = £a, f;^ = 0, 

R^j = ‘^515] + e-2" [rz/' - (l - ] a,-, 

Ck = {1- A) z>2 - 2ufj,'e^-^ + 

+A -z>e^-" - (2r/i' + 1) 

3 

Cv = J2^v\ (4-3) 

s=0 

where a prime denotes the partial derivative with respect to r, = 6 ^6j +sin^ 
and Cy^^s are given by Eq. (3.4). The Hamiltonian constraint (2.15) reads 

j (/li^ + £y - SttGJ*) eVdr = 0, (4.4) 

while the momentum constraint (2.17) yields 

(1 — A)|e^“'^ + 2(/rV — 1)] — z>'r^| + 2r — z>) 

= -STrGr^e-^+a, (4.5) 

where 

ji = e-(^+^)(n,0,0). 
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It can also be shown that Eqs. (2.19) and (2.20) now read 


(A^r^ - 1) + l] - 2(^z/' - 

+ (1 - A)|e2'^ (-r^h" + r^hV' - 2rz>') 

+e^+" [r2(/' + 3/iV" - - 3/i'V' + /i'^ - Sz/'/x'^ 

+2/i'z/'2) + 2r(2/i" - z/" + 2(z/' - /x')^) + 2z/'] } 

= SnCr^e^'^J^, (4.6) 

2rz/' - (A^r^ - l) + ij = ATiGr^e^'^J a- (4.7) 

The dynamical equations (2.22), on the other hand, yield 

(1 - A)r + fi' - z>') - + iz>^) + - /iV') 

+ 2 (/i' + Az/') + (4A - 3) ^ - 2e"+^ (A/1 + z>) 

+ = -r (^Frr + ^4 + j, (4.8) 

Ar^/i" - /i'z/'j + (2A - l)(2/i' - z/') + ^(3A + l)r/i'^ 

+ ^Az> + ^(A + l)z>^jre^'^ - [(2A - l)/i + r/i'(h + A/i) + Ar(z>' + /i')] 

2u 

= -— {Fee + + STTGr^p,) , (4.9) 

where 


Tij = e^''pr5l5] + r'^pe^ij, 

= ^(A' + Az/')(5[5} + e-2"p(A"-z/'A')+r(A' + Az/') 

and Fij is given by 


p2!/ 2 

1 — 
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o„— 2!^ p 

= - — [6e2" + - 8r^p" + 12r^ (lyf - tJ d[S; 

+ —^ + 4z/(3)r3 + 2Ar^ {u'f 

-2rv' + 14rV" + 7) - t] %, 

p-2v p 1 

(Fs)^,- =-^ [4e2" + - 6rV" + Qr^ (z/f - sj 515] 

0—4!^ r 

H-^ + 3z/^^V^ + 18r^ {v'f 

—rv' (—46^*^ + 21r^z/" + lO) — 5 VLij, 

[F^).. = - — (e2^ + 2rv' - l) ^226^^ + - 2Ar‘^v'' 

+40r2 {v'f - 2 (e^" - l) rv' - 23j 515] 

+ ^-^\ 240r^ (z/')^ + 4 (iSe^'^ - 17) {v'f 
-12r2 {v'f (-lle^^ + 22rV" + 15) 

+3rz/' - + le^'^ + 8z/®r^ - 28 - l) rV" - 1 

+2 12r^ {v''f - 24 - l) rV" 

+ {e^'' - l) (226^'^ + - 23) |%, 

9„—4;/ r 

(F 5 )ij =-^ eOr^ {v'f + {e^'^ - l) (iGe^^^ + - lAr'^v” - 17) 

+ (216^'^ - 17) {v'f - 4rz/' {-7e^'' + 9rV" + 7) 

9„—6j^ f 

+ 7^ J 18r^ (z/")^ + 180r^ {v'f + - l) (326^^ + 26^^^ + - 34) 

+21 {2e^'' - 1) {v'f - 28 - l) rV" - (z/')^ (-776^^^ + 198rV" + lOl) 

+rz/' 3 {—8e^'' + + 6z/*'^V^ + 3) — (496^^^ — 4l) r^v” 
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p-4!^ ( 

(FG)ij = \ -SOr^ {v'f - - l) {l3e^^ + - 6rV" - 14) 

{uf + 6rz/' {-2e^’' + br^u” + 2) |(5[(5; 

p-si' r 

+ 15r^ (z/")^ + 150r^ (z/')^ + - l) (26e^'^ + 26^“^ + Sz/^^V^ - 28) 

+ (iSe^'^ + 25) (z/')^ - 12 {e^'' - l) rV" - 3r^ (z/')^ (-lle^'^ + 55rV" + 12) 

+3rz/' - 126^^ + + 5z/®r=^ - - l) r^z/" + 8 

Op — iv ( 

{Fj)ij = —-—< — 2z/*'^V"^ + 15r"^ W)^ + 40r^ (z/')^ + 4rz/' (2e^'^ + 5z/^^V^ — 6) 

-2r2 {v'f [-3e^'' + 41rV" + 15) - 4 {e^'' - 3) rV" - 7|(5[5} 
o„-6!/ r 

-— <^ 126^^ + 26^^^ + z/®r® + 120r® iv'f + - Gz/^r^ 

+rz/' 246^^ + - IGz/^^V^ + 127r^ {u"f - 2 - 33) rV" - 57 

-rV" (86^'^ + 25z/(^V^ - 24) + (z/')^ (226^^^ + 101z/®r^ - 102) 

-2r‘^ {v'f (-6e^‘' + lG3rV" + 45) - 

p—4i' f 

(F8)ij = —G (2e^^ + ~ 3) + 45r^ (^")^ + 120r^ (ly')^ + lOr^ (z/')^ 

-8 - 4) rV" - (z/')^ (-12e^'^ + 24GrV" + 77) 

+2rz/' (86 ^^ + 30z/( 3V3 - 3rV" - 32) |(5[5)' + ^|-24e2'^ - 126^^^ - 3z/®r^ 

—3G0r^ (z/')^ — 3r^ W)^ — 30r^ {v')‘^ — 4e^'^z/*'^V^ + lGz/*^^V^ 

+rV" (iGe^'^ + 75 z/(3V^ - G4) + 2r^ {v'f {-12e^'' + 489rV" + 113) 

-r^ {Ff (446^'^ + 303z/(3V3 - 33r'^F' - 2G9) - rz/' 381r^ {F'f - 2 [Ue^'^ - 85) r^z/" 

+G (86^'^ + - 8z/(^V^ - 25) + 3p^^^r^ + 3G|ai. (4.10) 

where v' = dvjdr and Qij = SfSj + sin^ 9SfSj. 

We define a fiuid with pr = pe as a. perfect fluid, which in general allows energy- 
flow along a radial direction, i.e., v does not not necessarily vanish [9G]. 
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The energy conservation law (2.28) now reads 

J dr e'" pn + (ph + ^Pr) z> + 4 (h - vfi) - 2 (^Ja + vJa^ = 0, (4.11) 

while the momentum conservation law (2.29) yields 

vp! - [v' - pi) -^{v - Pr +Pe) - - e'"~^ v + v{2v - p) = 0(4.12) 

To relate the quantities J*, J* and to the ones often used in General Relativity, 
in addition to the normal vector dehned in Eq. (2.32), we also introduce the 
spacelike unit vectors and 0^ by 

0^ = rsin05j. (4.13) 

In terms of these four unit vectors, the energy-momentum tensor for an anisotropic 
fluid can be written as 

= Pun^rii, + q{nf,Xu + n^Xfi) + PrXt^Xu + Pe {0^,9^ + 4’^l4’u) , (4.14) 

where pu, q, Pr and pe denote, respectively, the energy density, the heat flow along 
the radial direction and the radial and tangential pressures, measured by the observer 
with the four-velocity n^. This decomposition is consistent with the quantities J* and 
J* dehned by 

pH = V = e^q. (4.15) 

It should be noted that the dehnitions of the energy density pn, the radial pressure 
Pr and the heat how q are diherent from the ones dehned in a comoving frame in 
General Relativity. 

4-2 Junction Conditions across the Surface of a Collapsing Sphere 

The surface S of a spherically symmetric collapsing star naturally divides the 

spacetime M into two regions, the internal and the external regions, denoted by M~ 
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Figure 4.1: The spacetime is divided into two regions, the internal M and external M"*", 
where M~ = ; r < and M+ = {x^ : r > The surface r = TZ{t) is denoted 

by S. 


and M~^ respectively, as shown schematically in Fig. 4.1. The surface S = dM = 
—dM^ is described by 

4>(t,r) = 0, (4.16) 

where <F(t,r) = r — TZ{t). The spherical symmetry implies that the ADM variables 
on M take the form (5.5). 

Preliminaries 

We assume that the normal vector to the hypersurface S with components 

= (4.17) 

+ (e^-^ + 77)5,^ 

is everywhere spacelike, i.e. 

This is the case if 77 is small enough. We may then define the vector field N = 
V$/||V$||g in a neighborhood of S. The vector field N should not be confused with 
the lapse function which in the present case is set to one, see Eq. (5.5). N has length 
one, i.e. NxN^ = 1, and the restriction of Al to S is the outward pointing unit normal 
vector field on S. 
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Let denote the Heaviside function defined by 


1 , $> 0 , 


//($) = 


- $ = 0 

0, <h < 0, 


(4.19) 


and let (5(<h) denote the delta distribution with support on E. By dehnition, (5(<h) 
acts on a smooth test function ip G C°°{M) of compact support by 


= (4.20) 

where dS = tArVol^ is the volume three-form induced by on S and ln denotes 
interior multiplication by N. The derivatives (5*^"^^(<h), n > 1, of d(*h) are dehned in a 
standard way and the following relations are valid [119]: 

dx^ ^ 

apr'S'”!*) = " = 0.1.2..... 

n = l,2,.... (4.21) 


If / is a function dehned in a neighborhood of E, we dehne the distribution 
J(5W(<|)) by letting it act on a test function (p by 

(/il”>(4),^) = (i‘”'(t)./4>). (4.22) 

The product /d(*h) is well dehned whenever / is and it depends only on the 
restriction /|s of / to E. More generally, the product is well dehned provided 

that / is C"' and it depends only on the values of / and its partial derivatives of order 
< n evaluated on E. 

Let F be a distribution on M of the form 

n 

F = F+//($) + F- [1 - F(<l>)] + F^”^(^)(5(^)(4>), (4.23) 

k=0 
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where the F^s are functions dehned in a neighborhood of S while F+ and F~ are 
sufficiently smooth functions defined on M+ and M~ respectively. We dehne the 
function F^ on M by 

+ F-[1(4.24) 
and we dehne the jump [F]“ of F across S by 

[F]~{x) = F^{x) — F~{x), a; e S. (4.25) 

We will also need the fact that the equation F = 0 is equivalent to the equations 

F^(a;) = 0, xeM^, (4.26) 


and 


~ ^)’ 4 ’ ^ plm{n-k) 

^ {j — k)\ 


k=0 


= 0 , 






where ^ acts on a function / by 

^ ^ yf /V 

||V4>||/‘ ’ 

and, more generally, for any j > 1, 

A proof of this fact is given at the end of the chapter. 
For n = 3, the conditions in (4.27) are 


plm(3)\^ = 0 , 


(9$ 

g2 plm(3) 

(9<F^ 

/q 3 plm{3) 

I ( 94>3 ■ 


= 0 , 


(9$ J 

q 2 plm{2) 


( 94>2 


+ 


Qplm(l) 

(9$ 


= 0 , 

s 

plm{0)'^ 


= 0 . 


(4.27) 


(4.28) 


(4.29) 


(4.30) 
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4-2.1.1. Proof of Eqs. (4.26)-(4.27). Let F be given by (4.23). We will show 
that the equation F = 0 is equivalent to the conditions (4.26) and (4.27). It is clear 
that the equation F = 0 is equivalent to (4.26) together with the condition 

n 

^ (<!>) = 0. (4.31) 

k=0 

It remains to show that (4.31) is equivalent to (4.27). 

Suppose hrst that (4.31) holds. Then, multiplying (4.31) by and using the 
recursion relation (4.21) repeatedly, we hnd 

G'5(^')($) = 0, 0<j<n, (4.32) 


where the function G is dehned in a neighborhood of S by 

n 

G{x) = (4.33) 

k=0 

Equation (4.32) with j = 0 implies that the restriction of G to S vanishes, i.e. 
G|s = 0. Equation (4.32) with j = I then gives 


0 = GF($) 


-4>F($) 

$ 




I.e. — 


0 . 


In terms of local coordinates {u^} such that = <h while the remaining coordinates 
{u^}j >2 parametrize the level surfaces of <h, we have 


0 = 


G 

¥ 


dG 


$=o 


Thus, G vanishes to the hrst order on S. Repeating the above procedure n times, we 
infer that G vanishes to the nth order on S: 


G 

<|)i 


= 0 i.e. 


d^G 




= 0, 0 < j < n. 


(4.34) 


< I>=0 
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The partial derivatives denoted in the local coordinates (u^) by ^ can be expressed 
invariantly as in (4.29). Substituting the expression (4.33) for G into (4.34), we hnd 


0 


n j 

fc=0 )'=0 


/j\ gj-rplm{k) gr^n-k 

\rj d¥~ 


i-l)^k\{n-k)\ 

k=n—j 


n 


gj-in-k) plm{k) 
gc^j-G-k) 


1 

^=0 


0 < j < n. 


Replacing k hy n — k, we hnd (4.27). 

Conversely, if (4.27) holds, then tracing the above steps backwards, we infer that 
(4.32), and hence also (4.31), holds. 

4-2.2 Distributional Metric Functions 

The held equations (2.15)-(2.22) involve second-order derivatives of the metric 
coefficients with respect to t and sixth-order derivatives with respect to x*. Thus, 
one might require that the metric coefficients be with respect to t and with 
respect to x*, where C"' indicates that the hrst n derivatives exist and are continuous 
across the hypersurface <l> = 0. However, this assumption eliminates the important 
case of an inhnitely thin shell of matter supported on S. Therefore, we will instead 
make weaker assumptions, so that a thin shell located on the hypersurface <h = 0 is in 
general allowed, and consider the case without a thin shell only as a particular case 
of our general treatment to be provided below. In fact, we shall impose the minimal 
requirement that the corresponding problem be mathematically meaningful in terms 
of distribution theory. Then, in review of Eqs. (4.4)-(4.12), we hnd that the cases 
A = 1 and A 7 ^ 1 have different dependencies on the derivatives of /i. In particular, 
the term p,'p!' appears when A 7 ^ 1. Thus, in the following we consider the two cases 
separately. 
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4-2.2.1. A = 1. In this case, we assume that: (a) fi and u are in each of the 
regions and M~ up to the boundary S; (b) /i is across S; (c) p is with 
respect to t and with respect to r across S. 

The above regularity assumptions ensure that the mathematically ill-dehned prod¬ 
ucts (5(<h)^ and (<h) do not appear in the held equations. Indeed, the terms in 

the held equations (4.4)-(4.12) that could lead to products of this type are 


, '2 
/i , 


fijj, 


V 


v"v"'. 


(4.35) 


Our assumptions imply that these terms may contain if (<h)^ but not (5(<h)^ or (5(<h)if (<h). 
In order to compute the derivatives of /i and z/, we note that 

H = fi+H {<!>)+{<!>)], 

u = = z/+if(4>)+ z/-[l-fi(4>)], (4.36) 


where the functions and are on M~^, while the functions fi~ and z/“ are 
on M~. Let Vs denote an open neighborhood of S. Let and z/"*" denote C^- 
extensions of jj,~^ and z/+ to U Vs. Let fl~ and z>“ denote O^-extensions of ii~ and 
z/“ to M~ U Vs. Then the functions 


fi = H — jj, 


V = V — V 


(4.37) 


are dehned on Vs and the following relations are valid on S whenever a -|- /i < 5: 


h = [h] , 


Qa+P 


Qa+ji - - 
Qa+P 

dt^dr^ 


(4.38) 

(4.39) 
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Since n is across S, we find 


/^,r (/^,»’) ) 

/i,rt = 

f^,rr ~ ^l^,rr) “ 1 “ 

H rrr ~ ^l^,rrr)^ 2/i j.j.(5($) -|- (*^)' 


(4.40) 


Since fi is across S, the derivatives of /r’*' and /i in any direction tangential to S 
must coincide when evaluated on S. In particular, since the vector U dehned by 

= 6^ + il5^, (4.41) 


is tangential to S (i.e. U^N\ = 0), we obtain 

= [h,t] + = 0, 


that is, 

fit — (4.42) 

after Eq. (5.30) is taken into account. Then, from Eq. (4.40) one hnds fi^tr = fJ-rt, as 
it is expected. 

Similarly, since u is across S, we also have 

0 = U^[ux\- = [ut]- + n[iy,r]~- ( 4 - 43 ) 

But [v^rY = 0, because v is assumed to be with respect to r. Thus = 0. 

Therefore, v is in fact across S. The same argument applied to z/j and z/,, now 
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implies that v is in fact across S. We find 


= 




= 

(^,r) , 


/ \D 

= 


z/(3) = 


yW = 

(;y(4))i5 + p(3)^($). 

= 

(;y(5))0^2z>(^)(5($)+Z>(3)y(4>), 


(4.44) 


where = d^u/dr'^. We emphasize that the expressions on the right-hand sides of 
(4.40) and (4.44) are independent of the extensions used to dehne fi and z> in (4.37), 
because the values of fi, z>, and their partial derivatives of order < 5 are uniquely 
prescribed on S in view of (4.39). 


4-2.2.2. The Junction Conditions. We will hnd the junction conditions across 
S by substituting the expressions (4.40) and (4.44) for the derivatives of fi and u into 
the held equations (4.4)-(4.12). 

Suppose that the energy density pn = —2J* has the form 

OD 

n=0 

where it is understood that only hnitely many of the are nonzero. Since, by 

(4.3), 


the Hamiltonian constraint (4.4) reads 


Cv — 


'r<R.(t) 


-|- Cy + AnGpfj) e'^r'^dr 


+ / i^K + + 47rGp^) e'^r'^dr 

J r>'JZ{t) 

oo 


(4.46) 


+4rf^(-l)' 


dn 


n=0 


dC 


r='R.{t) 


{P 


Im{n) u^2\ _ 


H 


e r 


‘)=0. 
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Table 4.1. A list of all field equations for A = 1. 


Variation 

w. r. t. 

Name of equation 

General 

version 

Spherically 
symmetric case 

Junction 

condition 

lapse N{t) 

Hamiltonian 

constraint 

(2,15) 

(4.4) 

(4.46) 

shift iV* 

Momentum 

constraint 

(2,17) 

(4.5) 

(4.47) 


- 

(2.19) 

(4.6) 

(4.47) 

gauge held A 

- 

(2.20) 

(4.7) 

(4.47) 

metric gij 

Dynamical 

equations 

(2.22) 

(4.8) and (4.9) 

(4.48) and (4.49) 

- 

Energy 

conservation law 

(2.28) 

(4.11) 

(4.50) 

- 

Momentum 
conservation law 

(2.29) 

(4.12) 

(4.51) 


The left-hand sides of Eqs. (4.5), (4.6) and (4.7) have no supports on the hyper¬ 
surface r = TZ{t). Thus, these equations remain unchanged in the regions M~^ and 
M“, while on the hypersurface S they yield 

V = v^, = Ja = {Ja)^. (4.47) 

In fact, in order to avoid that the ill-dehned product i7(<h)(5(<h) arises from the term 
Ja^' in (4.12), we will assume that Ja is . 

The gauge held A has dimension [A] = 4, so the action cannot contain terms like 
with n > 2, that is, it must be linear in A. We therefore assume that A has the 
form 

OO 

A(t,r) = A^ + 

n=0 
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It follows that 


A^r = (Ar-)^+ + + 

n=l 

A,rr = (Ar-r)^ + [2i,, + + [i + 

OD 

+ ^ + A^™(«-2)] 5W($). 

n=2 

Thus, 

= ^1 (A,rf + + [i + (5($) 

CXD 

+5: 

n=l 

Feh = (F4)“ + 5^F«t'’”'”’'5'”’(4'). 

n=0 

where 

= e~‘^'' [r‘^{A^rr)^ - u^rr‘^{A^r)^ + r{A^r)^ + rv^rA^ - (1 - , 

^A,/m(0) _ g-2. ^ 

+r(i + - (1 - , 

= e-2^ [r^{A + 2Af;"(°) + - r2z/,^(A^"^(o) + 

+r(yl^"^(o) + - (1 - , 

^A,rm(n) ^ ^-2,. ^^2 ('^/m(n-2) ^ - rV,, 

^^^^/m(n-l) ^ ^ - (1 - e^^), U > 2. 

From Eq. (4.10) we hnd that the functions contain no delta functions 

1 fip-4!^ 

(n)„. = 

{Fs)rr = 

{Fr)ee = {Fj)^ - Sre"®" [(2z>(^) - 16uM^'>)6{^) + z>(3)F($)], 

{Fs)e9 = {Fs)Z - 3re-6" [(2z>(^) - 16z/,,z>(3))5($) + z>(3)F($)]. 
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whereas 



Thus, (2.23) gives 


-Av 

Frr = {Frrf-{16gr + 69s)^i>^^^m, 

rpp 61 ^ 

Fee = -(8(77 +3(78)^[(2z>(')-16z/,,z>(3))5($) + z>(3)5'($)]. 


Writing pr in the form 


Vrit^r) = p^ + 

n=0 


we hnd that Eq. (4.8) remains unchanged in the regions and M , while on the 
hypersurface S it yields 


00 


E 

n=0 




— [e-2"(l - 2ru') + Agr^ - 


+r 




x5(”)($) = 0. (4.48) 

Using (4.27), Eq. (4.48) can be rewritten as a hierarchy of scalar equations on S. 

Similarly, Eq. (4.9) remains unchanged in the regions M+ and M“, while on the 
hypersurface S it yields 


+ 5^1 - 2rz/') + Agr^ - 

5W($) = 0. 


n=0 

^2v 


e 

H- 

r L 


plj(n) p pAIMn) p 


ee 


(4.49) 


Note that 

OD 

PH, = {PH,P + [p'A“> - np„]s{<s>) + - Kpr"-'h'”'(4). 

n=l 


and, by (4.47), 


v,t = (ut)'^ - 77h(5(<h). 
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Thus, in view of (4.47), the energy conservation law (4.11) takes the form 


that is. 


dre‘'r^i + [pI”,*”' “ 'Rph]'5(<*) 


n=l 

/ oo 


+4(p,)° + 4 E (4)) p, 

n=0 ' 

+4((nt)^-77h<5($)-n>t) 

-2((JA,i)^ + z/,t(J^)^)} =0, 


/ + / + «^,t(PH+ 4p,.) 

\Jr<n{t) Jr>n{t)/ 

+4nt - 4n/i,t - 2(7^,* + u^tJA))dr + 

-77pr^ + - 477h)l 

/ J r=Tl{t) 

+ E(-i)”|:;r eVfpSTP - Kp^”^' 


/m(n—1) 
H 


\r=n{t) ' 


+>',.(p7'”’+4p;’”'"'))] = 0 . 


(4.50) 


The momentum conservation law (4.12) remains unchanged in and M while 
on the hypersurface S it yields 

CXD 

-h5($) +p^5($) + ^[(p^”^(”)),^(5(”)($) 

n=0 

9 1 r 

+- -pj”^^"^)<5(")(4>) - -Ja {a + 7lf™W)<5(4>) 


^ J«(n) + ^(n) e^-'^77h(5(4>) = 0, 


where we have used that 
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This completes the general description of the junction conditions for the case A = 1, 
which are summarized in Table 1. 


4-2.2.3. A 7 ^ 1. In this case, the nonlinear terms 


,2 . / 
/i , nfi 


fj, p , 


z/ 


,w 2 


Z/'V'", 


(4.51) 


appear in the held equations (4.4)-(4.12). Thus, to ensure these held equations are 
well-dehned, we assume that: (a) /i and z/ are in each of the regions and M~ 
up to the boundary S; (b) jj, is with respect to t and with respect to r across 
S; (c) z/ is with respect to t and with respect to r across S. 

The same argument as above shows that u is and that p is across S. 
Equations (4.40) and (4.44) for the derivatives of p and v are still valid, but since fi 
now is C^, we have = /i,r = 0. It follows that all the junction conditions (4.46)~ 
(4.51) remain unchanged, except that the presence of the term /i'" in (4.6) implies 
that the expression for J^p now may include a delta function: 

+ (1 - A)—/.,„5(<1>). (4.52) 

In what follows, we will consider some specihc models of gravitational collapse 
for which the spacetime inside the collapsing sphere is described by the Friedman- 
Lemaitre-Robertson-Walker (FLRW) universe. 


4-3 Gravitational Collapse of Homogeneous and Isotropic Perfect Fluid 
In this section, we consider the gravitational collapse of a spherical cloud consist¬ 
ing of a homogeneous and isotropic perfect fluid described by the FLRW universe. 
Gravitational collapse of a homogeneous and isotropic dust fluid hlled in the whole 
space-time was considered in [120], using a method proposed in [121]: 

= —dP + — —r^ + r‘^d^Il\ , 

\ 1 — kr^ J 
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where /c = 0, ±1. Letting r = a{t)r, t = t, the corresponding ADM variables take 
the form (5.5) with N~ = 1, and 




jj, (t,r) = In 


-a{t)r 


.\/a? if) — 

where a < 0 for a collapsing cloud. For a perfect fluid, we assume that 


(4.53) 


Ve =Pr =P = 0- 


(4.54) 


We anticipate that the junction condition for u requires k = 0. Then, we hnd that 
u~{t,r) = 0, /i“(t,r) = In (—rif), (A; = 0), (4.55) 

where H = a{t)/a(t), and that 

£- = 3(1 - 3A)iL^ Cy = 2A, 

C~ = Cl = 0, = 2 A 34 I-. (4.56) 

It is easy to verify that the momentum constraint (4.5) is satisfied, whereas the 
equations (4.6) and (4.7) obtained by variation with respect to 99 and A respectively, 
reduce to 

?>kgH + S'KGJ- = 0, (4.57) 

AnGJ^ + ^9 = 0 - 

Since z/“ = 0, we have and the hrst dynamical equation (4.8) reduces 

to the condition 
4 

-a‘^A~ + 2a?AgA~ + 2(3A — l)ad + (3A — l)a^ + 2a^(87rGp“ — A) = 0. 

If this condition is satished the second dynamical equation (4.8) also holds provided 
that A~ — rA~j. = 0. On the other hand, the momentum conservation law (4.12) 
reduces to JA^~r = 0- We conclude that the general solution when A; = 0 is given by 

_ SAgH _ Ag 

^ 87rG ’ ^ 47rG’ 
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(4.58) 



with A = A (t) being given by 


AgA + (3-^ - 1) ( - + ^ 


— A = —SttGp 


(4.59) 


In the rest of this section, we consider only the case where Ag = 0. Then, Eq. (4.58) 
yields 


Ja = = 0, (4-60) 

for which Eq. (4.59) shows that now A~{t) is an arbitrary function of f, and a{t) is 
given by 

(3A — 1) H—— A = —87rGp~. (4.61) 

It is interesting to note that, since the Hamiltonian constraint is global, there is 
no analog of the Friedman equation in the current situation. This is in contrast to 
the case of HL cosmology [76], where a Friedman-like equation still exists, because 
of the homogeneity and isotropy of the whole universe. Although there is no analog 
of the Birkhoff theorem in HL theory, so that the spacetime outside the collapsing 
cloud can be either static or dynamical, we assume in this chapter that the exterior 
solution is a static spherically symmetric vacuum spacetime. We also assume that 
the value of A^ is the same in the exterior and interior regions, i.e. 

A+ = A; = 0. (4.62) 

It is convenient to consider the cases A = 1 and A 7 ^ 1 separately. 


4-3.1 Gravitational Collapse with A = 1 

We hrst consider the case of A = 1. In this case, the static spherically symmetric 
exterior vacuum solution has the form [78] 


= fi+{r) = iln 

Z/+ = 0, 


9rn+ 1 9 

-1—Ar^ — 2A'^(r) H— / A'^{r')dr' ), 

r 3 r 
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(4.63) 



for which we hnd that 


C+ = - 2A, £+ = 2A, £+ = 0, 

v+ = j+ = j+=p+= 0, (4.64) 

where m'^, tq are constants and = A'^{r) is a function of r only, yet to be 
determined. 

As mentioned previously, the condition that u be continuous across S implies that 
k = 0. We let the interior solution be of the form (4.55), and assume that the thin 
shell of matter separating the interior and exterior solutions is such that 

p = pA n = 0, = 4"^(°)5($), 

Pe=Pe + pH = Ph + pJr^°^(5(4>), 

A = A^ + A^”^(°)(5(<h), Ja = 0, (4.65) 


where Ph = = pt = pt = 0. 


Proposition 4.1. For the spacetime defined by (4-55), (4-61), (4-63), the six junction 
conditions (4-46)-(4-51) reduce to the following six conditions: 


6+ 2A + An +4 f A+rdr + 47rGpj^*'°V^ 

3 .JHit) 

= 0 , 


r=Tl(t) 


= 0,(4.66) 
(4.67) 


A(t,r) is continuous across S, 

Aj=n(^- hA ill - HU) - 87rGpl'^^^'> HU, 


(4.68) 

(4.69) 


pj"(°^(t,77(t)) pS"(°)(0,77(0))+ I(]H{s)n{sfpfi^fis) 


(4.70) 

(4.71) 




rp + 2pg™^*'°^ = 0 on S. 


Moreover, the condition that p he continuous across S implies that 




A-377^ 


nn-HHfTZ^ 


(4.72) 
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Proof. For the spacetime defined by (4.64)-(4.65), condition (4.46) reduces to 


{—QH{tY + 2K + AnGp^{t)) / + 4 / Afrdr + = 0, 

Jo Jn{t) ’ r='R.{t) 

which yields (4.66). Moreover, condition (4.47) reduces immediately to (4.67). 

Conditions (4.48) and (4.49) reduce to 




(4.73) 


r{p,^re^^ + 77/i,re^)d(<h) 

^A7m(2)y,(^) ^ = 0, 


(4.74) 


respectively, where we have used that 


Frr = {Fr, 


Fee = {Fei 


(4.75) 


F;t = -|(Ar)^+[^ + dlfr^°^]d(4>)+7l^”^®d'(4>)|, 
Fe'e = + 


(4.76) 


where 


= r‘^{A^rr)^ + r{A^r)^, 


^A,lm(0) 




r\2A^r + 7lf,A°^) + r(i + 

A(i + 27lfA°)) + 


?A,Im{2) _ ^2 j^lm(0) 


Thus, equation (4.73) can be written as 


[i + AfA°)](5($) + A^"^(°A'(4>) = 0. 


( 4 . 77 ) 



Thus, by (4.30), = 0. Hence, = 0 which gives 




Equation (4.77) then gives A|s = 0 so that in fact A is continuous across S, which 
proves (4.68). Equation (4.74) can now be written as 

/i,,(e2^ + Tte'^) + 2A^r + SttG'pJ’”® 1 (5($) + iy(4>) = 0. 

In view of (4.30) this yields 


/m(0) _ dA 


/i,,(e"'^ + 77e^) + 2H,, + SttCp^ = — on S. 


(4.78) 


Now observe that if a function f{t,r) is across S, then 


(9$ dr 


on S. 


(4.79) 


Indeed, the continuity of / implies that the derivative of / in any direction tangential 
to E must vanish when evaluated on S; thus f^t + T^f,r = 0 on S. A computation 
using (4.17), (5.142), and (4.28) now gives (4.79). 

On the other hand, since 




fr , 

’ r 


we hnd 


p,, = i(Ar-2At)e-2/^^--. 

2 ’ r 


(4.80) 


Inserting the equations (4.79) and (4.80) into (4.78), we hnd 


- (Ar - 2A+)e-2/^ - - ) (e^^^ + TZeA + A , 
2 ’ r / 


+87rGpg™'*'°^ = 0 on S. 


Since A^r = A+ = A^TZ simplihcation yields (4.69). 
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Condition (4.50) reduces to 


rn{t) 


e^r'^Pn^t^r + 


/m(0) ^ ^ 

PH,t - T^Ph 


A 

r=Tlit) dr 


r='R.{t) 




= 0 . 


That is, 




r-^dr + n{tf (t, 7^(t)) + 7^(t)p^(t) 


a{t) 


+27^(t)7^(t)pJ”®(^,7^(t)) +7^(t)27^(t)pJ”®(t,7^(t)) = 0. 


Consequently, 


4 dt 


+ 2^p;‘"’'(t,1Z(t)) + K(i)p„(«) = 0. 


TZ{t) lm{0). 


nit) 


(4.81) 


Solving this differential equation for we find (4.71). 

Condition (4.51) reduces to 


p + (5($) = 0. 


This yields (4.71). 

Finally, the condition that p be continuous across S can be written as 


+lA7^(t)2_2yl+(7^(^))+ ^ 


pTlp) 


A+(r')dr' = (4.82) 


nit) 3 ^ 77(t)^ 

Since A is continuous across S, we have A+(77(t)) = A~it). Hence, multiplying (4.82) 
by n and then differentiating with respect to t, we find 


hn^n - 2A-n = 2HH^tn-^ + sH^n^n. 


Solving this equation for A^, we find (4.72). 


The conditions (4.69) and (4.72) imply that 


- - hA in - Hn) - 87 rGpJ ”^®77 = —^77 - m,t 77 , 

2 J 2 
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□ 


(4.83) 



Hh + {2H^ + 2H^t - A)7^ - leTrGpJ”^® 


0 . 


Solving this equation for TZ{t) we find the following equation which expresses TZ{t) in 
terms of H{t) and the pressure on the shell; 

7^(t) = e-^oh-)'=^J7^(0) + 167^G' / ^ (4.84) 

I Jo H{s) } 

where I(t) is dehned by 

I = 2H+^-^. (4.85) 

4-3.2 Dust Collapse with A = 1 

Suppose now that the perfect fluid in the interior region consists of dust, i.e. 

Pj =Pj = 0- (4.86) 

Then, the condition (4.71) implies that 

= 0. (4.87) 


Solving equation (4.61) for a(t) we hud 


a{t) 


j oocoshs |^^(t-to)^, A 7^0, 

|^ao(to -1)2/3, A = 0, 


(4.88) 


where oq and to are constants. In the following, let us consider the cases A 7 ^ 0 and 
A = 0, separately. 


4 .3.2.1. A > 0. In this case, substituting the expression for a{t) into (4.85) we 


obtain 



and then (4.84) yields 


7^(t) = 7^o coshs - t)^ , (4.89) 

where TZq is a constant. Condition (4.72) now implies that A~f = 0, i.e. A~{t) = Aq 
for some constant Aq. Then, by (4.68), A'^{TZ{t)) = Aq. That is, A+(r) = Aq for all 
r such that r = TZ{t) for some t. Hence, the form of (4.89) implies that = Aq for 
all {t, r) in the exterior region. This gives 


A{t,r) = Ao. 


(4.90) 


Condition (4.66) now implies 

= —77-0 


[-6H^ + 2A + 4nGpH{t))n{t) 
127rG 

A + 7rG[l + cosh(\/^(to — t))]p]j{t) 


(4.91) 


67rGcosh3(^^^(to — t)) 

Substituting this into condition (4.71), or its equivalent form (4.81), we infer that 
p^it) satishes: 

TZq 


12 


cosh 3 




—77o\/^sinh 




(^0 - t)^ |piy,i(t) = 0, 


i.e. 


Pnit) = Ph\ 


(4.92) 


where is a constant. All the conditions of Proposition 4.1 are now satished. It 
only remains to consider the condition that p be continuous across S. This condition 
reduces to 

= ^ + iA772-2Ao + |Ao(77-ro)- 77^772 
6m+ — 6Aoro + TT-gA 
377o cosh^/^(^(to - 1 )) 

69 


0 



That is, the parameter tq is hxed by 


ro = 


6m+ + TZlA 
6 A() 


(4.93) 


This implies that 


, 1, /Ar2 A7^3 

" = 2 '"— 


(4.94) 


Since all the held equations and junction conditions are now satished we have proved 
the following result. 


Proposition 4.2. Hofava-Lifshitz gravity admits the following explicit solution when 
A = 1 and A > 0.' 


^ =ln(-hf(t)r), 

z/ = 0, H{t) =-^tanh(^^^^{to-t)^, 

n{t) = 7^o coshi - t)^ , 

Pr=Pe = 0, pH{t)=P%\ A(t,r) = Ao, 

is given by (4-91), 


(4.95) 


where to, TZq, Aq, and p® are constants and M+ = — A7^g/6. 


For t < to the dust cloud is contracting. As t —)■ to, the radius of the dust sphere 
approaches its minimal value of 7^ = TZq at t = to, and the function approaches 
zero: 

lim TZ{t) = TZo, lim = g, 

as shown schematically in Fig. 4.2. After the star collapses to this point, it is not 
clear how spacetime evolutes, because /i+ becomes unbounded as one can see from 
Eq. (4.95), for which the extrinsic scalar , 

+ -] , (4.96) 

r / 


K^{ 


ri = 


At(r) 
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Figure 4.2; The evolution of the surface of the collapsing star for A = 1 and A > 0, given 
by Eq. (4.95). At the moment t = to, the star collapses to its minimal radius TZ{to) = TZq, 
at which the extrinsic curvature becomes unbounded, while the four-dimensional Ricci 
scalar remains finite. 


also becomes unbounded, which indicates the existence of a scalar singularity at this 
point [122]. However, such a singularity is weak. In particular, the corresponding 
four-dimensional Ricci scalar remains hnite, = 4A. Thus, it is not clear whether 
the spacetime across this point is extendible or not. 

In addition, Eq. (4.91) shows that and pjj cannot both be positive. To 

understand this, letting M = —A7^q/ 6 we can write in the form 


1 , /2M Ar2 


= X In 

^ 2 \ r 


+ 


(4.97) 


However, this is nothing but the Schwarzschild-de Sitter solution with mass M and 
a cosmological constant A, where M is negative. 


4-3.2.2. A < 0. In this case, substituting the expression for a(t) into (4.85) we 
obtain 




and then (4.84) yields 


n{t) = TZo COS3 - to) ), 
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(4.98) 





where TZq is another constant. Condition (4.72) now implies that A~^ = 0, i.e. A~{t) = 
Aq for some constant Aq. Then, by (4.68), A'^{TZ{t)) = Aq. That is, A+(r) = Aq for 
all r such that r = TZ{t) for some t. We will assume that A'^ = Aq for all {t, r) in the 
exterior region, i.e. A{t,r) = Aq. Condition (4.66) now implies 

— (“677^ + 2A + 47rGpfj{t))7l{t) 


= 77-0 


127rG 

|A| - 7rG[l + cos(v/3|A|(t - to))]p7^(t) 


QnG cos 31 


-l^(t-to)) 


(4.99) 


Substituting this into condition (4.71), or its equivalent form (4.81), we infer that 
p'^it) satishes 

Pni^) ~ P^H 1 (4.100) 

where p® is a constant. All the conditions of Proposition 4.1 are now satished, while 
the condition that p be continuous across S reduces to 


0 = 


2imA 


+ AA77^ - 2Ao + ^Ao{n - tq) - H‘^n^ 


77 3 "77 

6m+ — 6Aoro + 77 qA 

377o cos2/^(^^^^^P^(t — to)) 


Thus, the parameter tq is hxed by 


ro = 


6m+ + 77|]A 
“6A^ 


This implies that 


^ 2 \ r 3 


(4.101) 


(4.102) 


where M = |A|77q/ 6. Clearly, this corresponds to the Schwarzschild-anti-de Sitter 

solution. For p’*' to be real, we must assume that r < 77o. Similar to the last case, the 

extrinsic curvature K~^ at r = 77o becomes unbounded, while the four-dimensional 

Ricci scalar *^^^7? remains constant. Thus, in this case it is also not clear whether or 

not the spacetime is extendible cross r = 77o. 

In any case, all the held equations and junction conditions are now satished for 

^ < 77o, and we have proved the following result. 
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Figure 4.3; The evolution of the surface of the collapsing star for A = 1 and A < 0, given 
by Eq. (4.103). The star starts to collapse at a time t = ti > to. At the later time t = tg, 
at which TZ{ts) = 0, the star collapses and a central singularity is formed. 


Proposition 4.3. Hofava-Lifshitz gravity admits the following explicit solution when 
A = 1 and A < 0.' 

p-=\n{-H{t)r), 

z/ = 0, H{t) =tan (^^l^^-{t-to)^, 

7^(t) = 7^ocosi - to)^, (4.103) 

Pr=Pe = 0, p-{t)=p%\ A{t,r)=Ao, 

is given by (4-99), 

where to, TZo, Ao, and p® are constants. 

The evolution of the surface of the collapsing star is illustrated in Fig. 4.3. The 
collapse starts at an initial time ti < to, and at time t = tg, the star collapses 
to a central singularity at which we have TZ{tg) = 0, where tg = to + t^/ a/ 3|A|. 
Equation (4.99) shows that now both and pjj can be positive, provided that 

|A| > 27rG'pg^ 
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4-3.2.3. A = 0. In this case, substituting the expression (4.88) for a(t) into 


(4.85) we obtain 


and then (4.84) yields 


I{t) 


2 

3(to -1)’ 


7^(^) =7^o(^o-t)^ (4.104) 

where TZq is a constant. Condition (4.72) now implies that = 0, i.e. A~(t) = Aq 
for some constant Aq. Then, by (4.68), A~^{TZ{t)) = Aq. That is, A+(r) = Aq for all 
r such that r = TZ{t) for some t. Hence (4.104) implies that A~^ = Aq for all {t,r) in 
the exterior region. Thus, in the present case we also have A{t,r) = Aq. Condition 
(4.66) now implies 

Ph IC 

^ „ 2-?>Gn{tQ-tfpH{t) 

° 9G7r(to — 

Substituting this into condition (4.71), or its equivalent form (4.81), we infer that 


PhY) — Ph\ 


(4.105) 


where is a constant. All the conditions of Proposition 4.1 are now satished, and 
the condition that p be continuous across E becomes 

0 = ^ - 2A„ + - r„) - 

^ ^ 9m+ - OApro - 2770 
977(t) 


Hence, the parameter tq is hxed to 


^^0 


9m+ — 277-0 

9l^ 


(4.106) 


which implies that 


A 




(4.107) 
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where Vg = 47^q/ 9. This is nothing but is the Schwarzschild solution written in the 
Painleve-Gullstrand coordinates [97]. All the held equations and junction conditions 
are satished, so we have proved the following result. 


Proposition 4.4. Hofava-Lifshitz gravity admits the following explicit solution when 
A = 1 and A = 0.- 


= fi = ln(-if(i)r), 

n{t) = 77o(to 

Pr=Pe = 0, pnif) = p%\ A{t, r) = Aq, 

Jmio) _ ^ 2 - 3G7r(to - tfpP 
~ ° 9G7r(to-t)A3 ’ 


(4.108) 


where to, TZq, Aq, and p^^ are constants. 


The evolution of the surface of the collapsing star is shown in Fig. 4.4. The star 
begins to collapse at the moment ti with a radius 77j[= until the moment t = to, 

at which we have TZ{to) = 0 and a central singularity is formed. The spacetime outside 
of the star is given by the Schwarzschild solution. Thus, as in GR, the Schwarzschild 
spacetime can be formed by the collapse of a homogeneous and isotropic dust perfect 
huid [108]. We note that > 0 for 


SGttp® 


<t < to 


4-3.3 Gravitational Collapse with A 7 ^ 1 

We now consider the case of A 7 ^ 1. For an exterior static spherically symmetric 
vacuum spacetime with A 7 ^ 1 and A^ = 0, equation (4.7) implies that 


u = 



(4.109) 
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Figure 4.4; The evolution of the surface of the collapsing star for A = 1 and A = 0, given 
by Eq. (4.108). At the moment t = ti < to, the star starts to collapse until the moment 
t = to, at which we have TZ{to) = 0, whereby a central singularity is formed. 


where B is a constant. On the other hand, for the interior FLRW region, we have 

‘'■=-^‘"( 1 -'=^)^ (4.110) 

Hence, the condition = z /7 on S implies that 0 = kTZ{t)‘^. Consequently, in order 
for a solution with TZ{t) 7 ^ 0 to exist, we must have k = 0. The conditions that u 
and be continuous across S then reduce to 2B/lZ{t) = 0. Thus, in order for a 
nontrivial solution to exist we must have k = B = 0. Thus, we have 


//“ = //■*■ = 0, fi~ = In (^—rUy (4.111) 

On the other hand, since A 7 ^ 1, the momentum constraint (4.5) yields 

/i+(r) = In ^Cir + , (4.112) 

where Ci and C 2 are constants. The held equations (4.6)-(4.9) are then satished 
provided that 


A+(r) 


^0 g^4 


+ 


3(1 - dA)^!^ + 2A 2 
8 ^ 


(4.113) 


where Aq is a constant. It is interesting to note that this class of solutions was hrst 

found in [81] in the IR limit. However, since the restriction of the spacetime to the 
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leaves t = constant is flat, we have Rij = 0 , and the higher-order derivative terms of 
Rij vanish identically, so they are also solutions of the full theory. Moreover, since 

+ _ -2C2 _ _ 1 

Cir^ + Car’ r’ 

we hnd 


, _ -3C'2 

Cir^ + C2r 

Thus, the requirement that /x is implies that C 2 
requires that H(t) = —Ci is a constant and so 


(4.114) 

0. The continuity of fi then 


a{t) = a^e 


It follows that is smooth across S. Note also that the asymptotical-flatness condition 
requires Ci = 0. However, in the following we leave the possibility of Ci 7 ^ 0 open. 
We hnd that 


£+= 3^2(1-3A), £+ = 2 A, £+ = 0, 

n+ = 0, J+ = 0, J+ = 0, p+=0. (4.115) 

In order for the integral over the exterior region in the Hamiltonian constraint (4.46) 
to converge, we also need to assume that 

3^12(1-3A) + 2A = 0. (4.116) 

Thus, A+(r) = Aq is a constant and equation (4.61) implies that p~{t) = 0, that is, 
the perfect huid in the interior region consists of dust. 

Similar to the case with A = 1, the interior solution is still of the form (4.60), i.e. 

= 0, A = A (t). 

In view of (4.56), we have 

= 0, £^ = 0, = 3(1 - 3A)R^, 

£y = 2A, V- = 0. (4.117) 
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We assume that the thin shell of matter separating the interior and exterior solutions 
is such that 


p^ = 0, v = 0, 

Pe = {<!>), = p- 

A = + Ja = 0, Ca = C-a, (4.118) 


with p^ = p^(t). 


Proposition 4.5. For the spacetime defined by (4-lll)-(4.118), the six junction con¬ 
ditions (4-46)-(4-51) reduce to the following six conditions: 


+p'A°’frn(t)) = 0 , 

j'™'”’ = 0 . 


A(t, r) = Aq is a constant, 

nt) Jm(0) 


= 0 on S 


d r 
dt 




+ 2 


7^(^) 


profit,n{t)) 


+7^(^)p^(^) +---= 0. 


(4.119) 

(4.120) 

(4.121) 

(4.122) 

(4.123) 


Proof. For the spacetime dehned by (4.111)-(4.118), condition (4.46) reduces to 


pn{t) 


(3(1 — 3A)C^ + 2A + 47rGpj:^(t)) J r^dr 
+ [ f3C'i2(l-3A)+ 2AV2dr 


JTlit) \ 


r=n{t) 


= 0 , 


which, in view of (4.116), yields (4.119). Moreover, equation (4.52) reduces to (4.120). 

The functions F and F^ are given by (4.75)-(4.76) also for A 7 ^ 1. Hence, condition 
(4.48) implies that A is continuous across S just like in the case of A = 1. Since 
A+ = A+ is constant and A [t) is independent of r, this gives (4.121). Condition 
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(4.49) then reduces to 


+ + 87rG'rpJ'"®5(<l>) = 0. 

Since A is a constant, this yields (4.122). 

Conditions (4.50) and (4.51) reduce to (4.123) and (4.122). □ 


Conditions (4.119) and (4.123) imply that 

fe(W)(Y^(«)-5) =0- 


Excluding the case of no collapse where Ti{t) is a constant, it follows that pjj must 
be a constant. If this is the case, then the junction conditions are satisfied provided 
that = —pJjTZ{t)/3. 

In summary, in the case A 7 ^ 1 a static spherical spacetime can be produced by 
gravitational collapse of a homogeneous and isotropic dust fluid. However, the space- 
time outside of such a fluid is not asymptotically flat, as one can see from Eqs. (4.111) 
and (4.112) with C 2 = 0. 


4-4 Conclusions 

In this chapter, we have studied gravitational collapse of a spherical cloud of fluid 

with a hnite radius in the framework of the nonrelativistic general covariant theory 

of HL gravity with the projectability condition and an arbitrary coupling constant A. 

Using distribution theory, we have developed the general junction conditions for such a 

collapsing spherical body, under the minimal requirement that the junctions should be 

mathematically meaningful in the sense of generalized functions. The general junction 

conditions have been summarized in Table I. 

As one of the simplest applications, we have studied a collapsing star that is made 

of a homogeneous and isotropic perfect fluid, while the external region is described 

by a stationary spacetime. We have found that the problem reduces to the matching 
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of six independent conditions that the Arnowitt-Deser-Misner variables {N, N’', Qij) 
and the gauge held A and Newtonian prepotential Lp must satisfy. 

For the case of a homogeneous and isotropic dust huid (a perfect huid with van¬ 
ishing pressure), we have found explicitly the space-time outside of the collapsing 
sphere. In particular, in the case A = 1, the external spacetimes are described by 
the Schwarzschild (anti-) de Sitter solutions, written in PainlevAGullstrand coordi¬ 
nates [97]. It is remarkable that the collapse of a homogeneous and isotropic dust to a 
Schwarzschild black hole, studied by Oppenheimer and Snyder in General Relativity 
more than 80 years ago [110], is a particular case. However, there are fundamental 
differences. First, in General Relativity a thin shell does not necessarily appear on 
the surface of the collapsing sphere [ 110 ], while in the current case we have shown 
that such a thin shell must exist, in General Relativity, because of the local conserva¬ 
tion of energy of the collapsing body, the energy density of the dust fluid is inversely 
proportional to the cube of the radius of the fluid, while in the current case it remains 
a constant, as the conservation law is global [cf. Eq. (2.15)] and the energy of the 
collapsing star is not necessarily conserved locally. 

In the case A 7 ^ 1, the space-time outside of the homogeneous and isotropic dust 
fluid is described by Eqs. (4.111) and (4.112) with C 2 = 0. It is clear that such a space- 
time is not asymptotically flat. Therefore, in this case to obtain an asymptotically 
flat space-time outside of a collapsing dust fluid, it must not be homogeneous and/or 
isotropic. 

From the above simple examples, one can already see the signihcant differences 
between the HL theory and General Relativity in the strong gravitational held regime. 
Therefore, it is very interesting to study gravitational collapse of more general huids, 
such as perfect huids with diherent equations of state, or anisotropic huids with or 
without heat hows. Particular attention should be paid to the roles that the equation 
of state and heat hows might play. It would be extremely interesting to study the 
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implications for black hole physics, or more generally for (observational) astrophysics 
and cosmology [108]. Since the general formulas have been already laid down in this 
chapter, we expect that such studies can be carried out easily. 

As emphasized previously, our treatment of the junction conditions of a collapsing 
star presented in this chapter can be easily generalized to other versions of Hofava- 
Lifshitz gravity or, more generally, to any model of a higher-order derivative gravity 
theory. 
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CHAPTER FIVE 


Global Structure of Spacetime 

This chapter published as: J. Greenwald, J. Lenells, J. X. Lu, V. H. Satheeshkumar 
and A. Wang, “Black holes and global structures of spherical spacetimes in 
Hofava-Lifshitz Theory,” Phys. Rev. D 84, 084040 (2011). 


5 .1 Introduction 

In the HL theory, due to the breaking of the general covariance, the dispersion 
relations of particles usually contain high order momentum terms [123, 124, 125, 126, 
127], 


z-l 


Ui. = m 


1 + ^ ^ 


n=l 


for which the group velocity is given by [128] 



(5.1) 



As an immediate result, the speed of light becomes unbounded in the UV. This makes 
the causal structure of the spacetimes quite different from that given in GR, where 
the light cone of a given point p plays a fundamental role in determining the causal 
relationship of p to other events [cf. Fig. 5.1]. However, once the general covariance 
is broken, the causal structure will be dramatically changed. For example, in the 
Newtonian theory, time is absolute and the speeds of signals are not limited. Then, 
the causal structure of a given point p is uniquely determined by the time difference. 
At = tp — tq, between the two events. In particular, if At > 0, the event q is to the 
past of p; if At < 0, it is to the future; and if At = 0, the two events are simultaneous. 

Another consequence of the breaking of the general covariance is that free particles 

now do not follow geodesics. This immediately makes all the definitions of black holes 

given in General Relativity invalid [129, 130, 131, 132]. To provide a proper definition 

of black holes, anisotropic conformal boundaries [133] and kinematics of particles 
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Figure 5.1; (a) The light cone of the event p in special relativity, (b) The causal structure 
of the point p in Newtonian theory. 


[134] have been studied within the HL framework. In this chapter, we shall adopt 
the approach of Kiritsis and Kohnas (KK) [135], where a horizon is dehned as the 
inhnitely redshifted 2-dimensional (closed) surface of massless test particles. Clearly, 
such a dehnition reduces to that given in General Relativity when the dispersion 
relation is relativistic [Where A„ = 0, as shown in Eq. (5.1).]. 

It should be noted that black holes in the HL theory with or without the pro- 
jectability condition have been extensively studied, mainly using the dehnition bor¬ 
rowed directly from GR. In this chapter, we shall show explicitly how these dehnitions 
are changed by considering some particular examples, found in the HMT set up with 
A = 1. 

Another interesting approach is the equivalence between the HL theory (without 
the projectability condition) and the Einstein-aether theory in the IR [136], where 
the former is equivalent to the latter for the case where the aether vector held is 
hyper surface-orthogonal. In the spherically symmetric case, this is not a restriction 
as the aether held now is always hypersurface-orthogonal. From such studies one 
already sees the difficulties to dehne black holes, because of the fact that diherent 
modes may have diherent velocities even in the IR. In [136], black holes are dehned 
to possess both a metric horizon and a spin-0 mode horizon. Since the equivalence 
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holds only in the IR, it is still unclear how to extend such dehnitions to high energy 
scales, where high order curvature terms become important. 


5.2 Black Holes in HL Theory 

Kristis and Kohnas (KK) considered a scalar field with a given dispersion relation 
F(C) [135]. In the geometrical optical approximations, ( is given by C = dijk^k^, 
where ki denotes the 3-momentum of the corresponding spin-0 particle. With this 
approximation, the trajectory of a test particle is given by 

_ 

Sp — / jTpdr 

Jo 



f -2 r 


t^ + e 

Jo 1 

[ ^ 


F(0 - 2(F'(0 


(5.3) 


where e is a one-dimensional einbein, and ( is now considered as a functional of 
t, x\ i, H and e, given by the relation. 


C |f’'(C)l" = + N‘i) (r + NH), (5.4) 

with t = dt/dr, etc. For detail, we refer readers to [135]. 

It should be noted that KK obtained the above action starting from a scalar 
field. So, strictly speaking, it is valid only for spin-0 test particles. However, what 
is really important in their derivations is the dispersion relationship F{(). As shown 
in [137], a spin-2 particle has a similar dispersion relation. It is expected that a spin-1 
test particle, such as photons, should have a similar dispersion relation too [128, 135]. 
Therefore, in the rest of this chapter and without proof, we simply consider the action 
(5.3) to describe all massless test particles. 

Spherically symmetric static spacetimes in the framework of the HMT setup were 
studied systematically in [79, 78], and the metric for static spherically symmetric 
spacetimes that preserve the form of Eq. (2.2) with the projectability condition can 
be cast in the form [103] given below. Note the slight difference between the gtr term 
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defined here and the one dehned in [103, 78]. 

ds'^ = —(?dt^ + [dr + e^~^cdtY + r'^d'^Q, (5.5) 

where d^Vt = d9‘^ + sin^ 9d(l)^, and 

jj, = jj,{r), v = v{r)^ TV* = 0,0} . (5.6) 

The corresponding timelike Killing vector is ^ = dt, and the diagonal case TV*" = 0 
corresponds to /i = —oo. 

However, to study black hole solutions in a more general case, in this (and only 
in this) section, we also consider the case without projectability condition, and write 
the metric as, 

ds^ = —N'^c^dt^ + j {dr + N^cdt)'^ + r‘^d^VL. (5.7) 

where iV, / and iV’’ are all functions of r. Without loss of generality, in the rest of 
the chapter we shall set c = 1, which is equivalent to the coordinate transformations 
Xq = ct, TV’’ = TV’’/c. Taking 


F{C)=C, (n = l,2,...). 


(5.8) 


Eq. (5.4) yields. 


c = 


r + N^t 


2/(2n-l) 


-) 

eV 


l/(2n-l) 


(5.9) 


neXJ 

Inserting this into Eq. (5.3), we hnd that, for radially moving particles, Cp is given 
by 

Cp = + ^(1 - 2n)eF^^-^F'pn/i2n-i)_ 

Then, from the equation 5Cp/5e = 0 we obtain 


^2^2 _ g2(n-l)/(2n-l)pn/(2n-l) _ 


(5.11) 


On the other hand, since SCp/St = 0, the Euler-Lagrange equation, 

^Fp _]_ / dCp \ _ ^ 

5t dr \ 5i ) ’ 
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yields 

NH - e2(-l)/(2n-l)i^pl/[2(2n-l)] _ ^ ^5_^2) 

V7 

where E is an integration constant, representing the total energy of the test particle. 

To solve Eqs. (5.11) and (5.12), we hrst consider the case n = 1, which corresponds 
to the relativistic dispersion relation. From such considerations, we shall see how to 
generalize the dehnition of black holes given in General Relativity to the HL theory 
where n is generically different from 1, as required by the renormalizability condition 
in the UV. 


5.2.1 n = l 

In this case, Eqs. (5.11) and (5.12) reduce, respectively, to. 


nH'^ - V 


nH - iv" 



0 , 

eE. 


(5.13) 

(5.14) 


Eq. (5.13) simply tells us that now the particle moves along null geodesics. The above 
equations can be easily solved according to whether iV’’ vanishes or not. 


5.2.I.I. = D. When 


0, from Eq. (5.13) we hnd 
dr 


dt = E 


iVx/7’ 


(5.15) 


where “+” (“—”) corresponds to out-going (in-going) light rays. If / has an a-th 
order zero and a b-th order zero at a surface, say, r = Vg, that is. 


f = fo{r){r -Tg^, N = No{r){r -Tgf^^, (5.16) 


where iVo(rg) ^ 0 and foi'f'g) 7^ 0, then from the above equations we hnd that in the 
neighborhood of r = r^. 


t ~ tn ± 



(r — VgY (“+^)/2^ a-I-6 7 ^ 2, 


No^o 


a + b = 2. 


(5.17) 



Therefore, when 


a + b>2,{n = l), (5.18) 

t becomes unbounded, as r ^ r^, at which the light rays are inhnitely redshifted. This 
indicates that an event horizon might exist at r = r^, provided that the spacetime has 
no curvature singularity there. A simple example is the Schwarzschild solution, = 
f = {r — rg)/r, which is also a solution of the HL theory without the projectability 
condition, but with the detailed balance condition softly broken [138], and for which 
we have a = b = 1. Clearly, it satishes the above condition with the equality, so 
r = Tg indeed dehnes a horizon. 


5.2.1.2. 7 ^ 0. When W 7 ^ 0, Eq. (5.13) yields 

f edr 

^ ^ ^ y iVyj-eiV-’ 

where e = +l(e = —1) corresponds to out-going (in-going) light rays. If 

H(r) = N./J-(N\ 


(5.19) 


(5.20) 


has 5-th order zero at Vg, 


H{r) = Ho{r){r - VgY, (5.21) 


with HQ{rg) 7 ^ 0, we hud that in the neighborhood r = Vg Eq. (5.19) yields 


t — to + 


Mr,) 


Tj))' - r,)* 5 1, 

ln(r — Tg), 5 = 1. 


(5.22) 


Clearly, when 


5 > 1, (n = 1), 


(5.23) 


|t| becomes unbounded as r ^ r^, and an event horizon might exist. 

The Schwarzschild solution in the Painleve-Gullstrand coordinates [97] is given 
by 

NlA = fs,k^l. = 


87 


(5.24) 



where ei = ±1. As shown in [79, 78], this is also a vacuum solution of the HL theory in 
the HMT setup [73]. Then, we find that H{r) = 1 — eie^rgjr. Thus, for the solution 
with €\ = +1, the time of the out-going null rays, measured by asymptotically flat 
observers, becomes unbounded at r^, and for the solution with e\ = —1, the time of 
the in-going null rays becomes unbounded. Therefore, an event horizon is indicated 
to exist at r = in both cases. 

In review of the above, KK generalized the notion of black holes dehned in General 
Relativity to the case of a non-standard dispersion relation [135]. In summary, a 
horizon is defined as a surface on which light rays are infinitely redshifted. It should 
be noted that this redshift should be understood as measured by asymptotically flat 
observers at A^(r 3> rg) ~ 1 and N^{r 3> rg) ~ 0, with r being the geometric radius, 
r = a/A/ 47r, of the 2-sphere: t,r = constants, where A denotes the area of the 
2-sphere. 


5.2.2 n>2 

In this case, eliminating e from Eqs. (5.11) and (5.12) we hnd that 


-p{r)X -q{r,E) = 0, 


(5.25) 


where 


X 

p(r) 


— \ l/(n—1) 

Vv\ _ Hr' + N'' 
i ) V ny/J 


77 ’ 


q{r,E) = 




(5.26) 


with r' = r/t = dr/dt. To solve the above equation, again it is found convenient to 
consider the cases = 0 and 0 separately. 


5.2.2.1. = 0. When = 0, Eq. (5.25) has the solution. 


dr 


t — to T c 




( 5 . 27 ) 



where e = +1 corresponds to outgoing rays, and e = —1 to ingoing rays. Thus, if / has 
an a-th order zero and a b-th order zero at r = Vg, as given by Eq. (5.16), we have 
(^r — Then, from the above, we hnd that the time t, measured 

by asymptotically flat observers, becomes inhnitely large at r = r^, provided that 
[135] 

b 

a + ->2. 5.28 

n 

For the solutions with the projectability condition [N = 1, 5 = 0), this is possible 
only when a >2. 

Considering again the Schwarzschild solution, = / = {'f' — rg)/r, one hnds that 
this does not satisfy the condition (5.28) with n > 2. Therefore, the Schwarzschild 
black hole in General Relativity is no longer a black hole in the HL theory, because 
of the non-relativistic dispersion relations (5.1). This is expected, since even in Gen¬ 
eral Relativity when quantum effects are taken into account, such as the Hawking 
radiation, classical black holes are no longer black. 

5.2.2.2. N'' 7 ^ 0. In this case, let us consider an ingoing ray r' < 0. Suppose 
there is a horizon located at r = th- Then r'(r) ~ 0 as we approach the horizon. 
Thus, if iV'’ > 0 and bounded away from zero, (r' -|- N^) will also be positive, when 
the ray is sufficiently near the horizon. Conversely, if < 0 and bounded away 
from zero, then (r' -|- N^) will be negative sufficiently near the horizon. Dehning H 
by H{r, E) = r', we hnd that for an ingoing ray near the horizon we have, 

« = <0 + /« fe . ( 5 . 29 ) 

H(r, E) = 6n77A'”-‘ - N', (5.30) 

where 

|l, > 0, 


- 1 , < 0 . 


(5.31) 



Dividing (5.25) by X and solving for X” we obtain 

Substituting this into (5.30), we find, 


H = {en- 1)X" + en^/f 


EN"-^ 


(5.32) 


It follows that if H has a zero at r = th, then 


eriy/jEN"- 
^ “ (en-l)X- 


(5.33) 


The expression on the rhs is positive (negative) for e = —1 (e = 1). Thus, H can have 
a zero only if e = —1. Thus, we will henceforth consider only this case. Differentiation 
of (5.32) with respect to r yields 

, ,— EX‘‘^~^ ^ h/ Af n —1 ,— Af n—l 

H'{r) = -in + m' - /' + ( 5 . 34 ) 

On the other hand, differentiation of (5.25) with respect to r yields 


_ 1 / 1 dN^ N'' df\ EN^-UN 

nX^-^ — ^ \i/7 dr 2 / 3/2 dr / ^ — i df 

Substituting the above into Eq. (5.34), we hnd that 


n+l(H, N^f , 2N^N' , 

2 [h 2 f ^ N-nN^ 


where 


(5.35) 


(5.36) 


Hi = 2E{n + l)fN^N’ + E{n - l)nN (^N^f - 2fN^'^ , 
H 2 ^ (n - l)nfN E+(n+ l)N^~ , 
If H has a zero of order d > 0 at we can write it in the form. 


(5.37) 


H{r) = Ho{rH){r - ThY H-, 


(5.38) 
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as r ^ th, where HqIth) 7 ^ 0. Therefore, 


H'ir] 


= < 


r=rH 


0 , (5>i, 

Hoirn), 6 = 1, 
±cxD, 0 < (5 < 1. 


Now t ^ 00 as r ^ ii and only if 


5 > 1 , 


(5.39) 


(5.40) 


which happens if and only if dH /is hnite. This gives an explicit condition on 
/, N, TV'’, E, n for the blow-up of t at th. 

It should be noted that th usually depends on the energy E of the test particles, 
as can be seen from the above and specihc examples considered below. 

Case n = 2: In this case, we have 


where 


H'{r) 


r=rH 


Hs 

2N[AEfN + 3{N^yy 


(5.41) 


H 3 = 3 4B1V= (iV 7' - 2/iV'') + SEfNN^N' - 3(N' fN' , (n = 2). (5.42) 


Again, for the Schwarzschild solution (5.24), we have 


X(r) = 


H{r) = 


+ 


y^AEr + ■ 


—3^rg{AEr + Vg) + AEr + Svg 


/rvg — ^Jr{AEr + Vg) 


so that H{r-) =0 has the solution. 


(5.43) 


Th = % (n = 2 ), 


(5.44) 


at which we have 


H\rH) = 


2 ^ 3/2 

\/3r„ 


(5.45) 
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Then, according to Eq. (5.39), we have 5 = 1, i.e., t diverges logarithmically as 
r ^ rjj. Therefore, in this case there does exist a horizon. But, the location of it 
depends on the energy E of the test particle, and approaches zero when E ^ Vg. This 
is understandable, as the speed of light is unbounded in the UV, and in principle the 
singularity located at r = 0 can be seen by asymptotically flat observers, as long as 
the light rays sent by the observers have sufficiently high energies. 

Case n = 3: In this case, we have 

— p(r)X — g(r, E) = 0. (5.46) 


Assuming that if(r) = 0 has a real and positive root th-, we hnd that 


H'{r) 


r=rH 


_ ih _ 

3N {27E^p/'^N - 16(iV’-)3) ’ 


where 


( 5 . 47 ) 


Hi = + 32{N'')‘‘N' - 162E^f/^N hNN'’ - /V/V') . (5.48) 

For the Schwarzschild solution (5.24), we have p{r) = —pTg/r and q{r,E) = E. 
Then, we hnd that 

+ y^X - E = 0, (5.49) 

from which we hnd that H{r) = 0 has a solution, 

/ \ 2/3 

= Tg , (n = 3), (5.51) 

which also depends on E, and approaches zero as E —)■ oo. Substituting rn into 
Eq. (5.47), we hnd H'^rn) = —27E'^/{16rg). That is, the hypersurface r = rn is also 
an observer-dependent horizon in the case n = 3, and the radius of the horizon is 
inversely proportional to the energy of the test particle. For E 3> r^, we have th — 0. 
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Another (simpler) consideration for the existence of the horizon is given as follows: 


First, from Eq. (5.26) we find that 


X = 


\^' _|_ ]\fr\ \ 




r \ 


eN^ \ 

ny/f) 


X 1 + 


H 


in - l)N^ J ’ 


(5.52) 


for r ~ r'H- Inserting it into Eq. (5.32), we have, to leading order. 


/ 


1 + 




\ 


\ ("- i )( Si ) 

Then, we obtain 


n/(n—1) 


/ 


H{r, E) = (en - l)iV" 


J 


1 + 








n/(n—1) 


15.53) 




{n + 1) 


-N^ 

«-vT 


n/(n—1) 


= 1 . 


(5.54) 


r=rH 


Given this, we can further simplify Eq. (5.53) to, 

/ 


2n 

n — 1 


H{r,E) = -{n + l)N^{rH) 


\ 


1 - 


EXi/("-i)(r) 


V 


(n + 1) 


-N^ir) \ 


n/(n—1) 


(5.55) 


J 


T'y/IiE J 

Then, using Eq. (5.38), we have the following constraint for N,N'',f to satisfy so 
that a horizon can indeed exist. 


EN^/i^-^)(r) 


{n + 1) 


\ n/{n-l) 
-N^ir) \ 


1 + 


2nifo(rj/) 

(n2 - l)N'-{rH) 


thY + • ■ • • 


(5.56) 


This equation can be first used to determine th and then 5, once X, N'^ and / are 
given. To illustrate how to use it, let us consider the Schwarzschild metric (5.24). For 
n = 2, th can be obtained simply from the above, and is given exactly by Eq. (5.44), 
for which we have 


that is, (5 = 1. 


EXi/("-i)(r) 


(n + 1) 


-Ny(r) 

«\/7W 


n/ (n—l) 


~ 1 + 


r — Th 
Th 


(5.57) 
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For n = 3, from Eq. (5.54) we find that th is given by Eq. (5.51), and 


EA^i/('^-i)(r) 


33/2 ^^ 3/4 


/ I 1 \ / —N'^(r) 


n/(n—1) 


4r. 


3/4 


~ 1 + -—(r - Th) H- 


4r 


(5.58) 


H 


Therefore, in this case we have <5 = 1 too. 

It should be noted that in the above analysis, we assumed that F{() = In more 
realistic models, the dispersion relation is a polynomial of C, as shown by Eq. (5.1), 
or more specihcally, 

where Ma and Mb are the energy scales, which can be signihcantly different from the 
Planck one [139]. Therefore, for observers in low energy scales, where C, Ma,Mb, 
the hrst term dominates, and some solutions, including the Schwarzschild solution, 
look like black holes, as shown in the case n = 1. But, for observers with high energies, 
those solutions may not be black holes any longer. Even if they are, their horizons in 
general are observer-dependent, as shown in the cases n = 2 and n = 3 explicitly for 
the Schwarzschild solution. To illustrate the main properties of the dispersion relation 
(5.59), we shall consider the case where only the hrst two terms are important. 


5.2.3 Trajectories of Test Particles with the Dispersion Relation F{Q = C + C^/^1 
For the sake of simplicity, we restrict ourselves to the case iV’’ = 0. Substituting 

f’(0 = C+^. (5.60) 

into Eq. (5.4), we hnd 

4 . 61 ) 

Solving this equation directly for (j yields a very complicated expression, and it is not 
clear how to proceed along this direction. Instead, we note that our goal is to hnd 
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the analog of equation (5.11), i.e. of the equation 5Cp/5e = 0, where 


r = 


1 r 

^[—i^ + e[F{0-2CF'{0]j 

2 iv' ^+Mf )• 


(5.62) 


Thus, we will hrst calculate S(/6e, implicitly by applying 5/5eto both sides of (5.61), 


which yields. 


2 Mjr 2 


& e3/(Ml + 8MlC+12C2)- 

Substituting this into the expression 

6e 2 V e2 Ml\ ^ [de Ml 5e 

we hnd the following analog of equation (5.11), 


(5.63) 


C {e^Ml + 2NH‘^) + 5e\2 + _ _ J = Q 


(5.64) 


where C, is given implicitly by Eq. (5.61). Note that in the limit Ma —t oo, the above 
equation reduces precisely to Eq. (5.11) for F{C,) = C, and iV^ = 0, as expected. 

On the other hand, the analog of Eq. (5.12) is simply 


NH = eE. 


(5.65) 


Using Eqs. (5.61) and (5.65) to eliminate r and t from Eq. (5.64), we find. 


& 2 Ml Af 2 ) 


In the limit Ma —t oo, this equation reduces to ( — ^ =0, which is again consistent 
with the case F{() = (. Solving this equation for (, we infer that 

Ml , Ma 


C = -^ + ^\/^e^ + mIn^. 


Substitution of this expression into Eq. (5.64) yields 


Ma [§J + Ml)+ 4 E^) ^2^2 3^2 

NJAE^ + MlN^ e 2 m 


-Ml=0. 


(5.66) 



Replacing e by N‘^t/E and then solving the resulting equation for r/t, we hud 

^2 fN {4:E^ + MlN‘^) (^y/IPTMp^ - MaN^ 

P ~ 2EmA ■ 

In the limit Ma —t oo, this equation becomes ^ = /iV^, which is again consistent 
with the case E{() = (. Thus, the trajectory is given by 

dr 


t — to + 


H{r, E) ’ 


(5.67) 


where 


H(t,E) = + Mlir- - M^N. ( 5 , 68 ) 

As an example, let us consider the Schwarzschild solution, = / = 1 — Vg/r, for 
which we hnd 

:{r - Vgf/^ + 0(^{r - rg)^/^y 


H = 2, 


MatT 


as r —)■ Tg, so that t remains hnite. On the other hand, as Ma —t oo, 




3^2 


M\ 


Thus, if we take the limit Ma —t cxd before letting the trajectory approach r^, then t 
will blow up logarithmically as r —)■ r^. As a result, a horizon exists in this limit. 
More generally, if / has an ath order zero and has a 6th order zero at r = r^. 


as given in Eq. (5.16), then, we hnd that 
H = 2 


lEfo{rg)No{r 


Ma 


_ rir, \ 2 + 4 O I (r — Ta) 2 + ^4 


as r —)■ Tg. It follows that 


t ~ to + 




+ 1 ^ 1 , 


EMrg)No{rg) 


Ma 


? 2 ' 4 

2 4 


ln(r - r,), f + | = 1. 


( 5 . 69 ) 


Therefore, t blows up as r ^ rg, if and only if 


®+5>2, 


which is exactly Eq. (5.28) for n = 2, as expected. 
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5.3 Vacuum Solutions with iV’’ = 0 


When = 0, the vacuum equations with J* = v = pr = pe = Ja = ^ yield 

the following most general solutions [78], 

/(r) = 1 + - - iV = l, N^ = 0 = p, (5.71) 

V o 

with the Hamiltonian constraint 

J Cve'^r^dr = 0, (5.72) 

where £y = Cv{r-,Kg,C,gs), as dehned in Eq. (2.13). 

The gauge held A must satisfy the equations, 


A' + An' -rFrr = 0, 

(5.73) 

r^{A” - n'A') +r{A’ + n’A) - A{l - 


+ e'^'^Fee = 0 , 

(5.74) 


where Fij is given by Eqs. (2.23) and (2.25). Then, from Eq. (5.73) we hud that 


A = Aoe-‘' 






(5.75) 


where Aq is an integration constant. The solutions with = 0 were hrst studied 
in [73, 79]. 

Since now we have iV = 1 and 5 = 0, Eq. (5.28) shows that a horizon exists 
only when a > 2. It can be shown that for the solutions given by Eq. (5.71), this 
is impossible for any chosen C and Kg. Therefore, it is concluded that the solutions 
given by Eq. (5.71) do not represent black holes. 

However, in some cases /(r) = 0 does have a real and positive root. So, there in¬ 
deed exists some kind of coordinate singularities, and to obtain a maximally (geodesi- 
cally) complete spacetime, some kind of extensions are needed. Because of the 
breaking of the general covariance and the restricted diffeomorphism ( 2 . 2 ), it is 

not clear if this requirement is still applicable here in the HL theory. Even if it 

97 



is not, some kind of extensions still seems needed. Such extensions are also needed 
in order to determine the range of r, from which the Hamiltonian constraint (5.72) 
can be carried out. Once this constraint is satished, one can integrate Eq. (5.75) 
to obtain the gauge held A. To this end, we divide the solutions into the cases: 
(i) C = Ag = 0, (n) C ^ 0, Ag = 0, (in) 0 = 0, 7 ^ 0, and (iv) O 7 ^ 0, 7 ^ 0. 

The hrst case is trivial, and it corresponds to the Minkowski spacetime with 1 / = A = 0 
and A = Aq. Thus, in the following we shall consider only the last three cases. 


5.3.1 0^0, A 3 = 0 

In this case the metric takes the form 

H ———c + ^^d^O, (5.76) 

r 

from which we hnd that 


Cv 


o. , 3^7302 3g,C^ 45gsC^ 



(5.77) 


where A = g^^C^ j2. To consider the Hamiltonian constraint (5.72), we need to further 
distinguish the cases O > 0 and O < 0. 


5.3.1.1. O > 0. When O > 0, the metric (5.76) is singular only at r = 0, so 
the solution covers the whole spacetime r G (0, cxd). The singularity at the center is 
a curvature one [ 122 ], as it can be seen from the expressions, 

302 




= 




2r6 ’ 

30=^ 

4^9 ’ 
4502 

2r8 


1 + - 
r 


(5.78) 


Since event horizons do not exist for O > 0, this singularity is also naked. Inserting 
it into Eq. (5.72), we hnd that the Hamiltonian constraint is satished only when 


A = 93 = Je = 9s = 0 . 


(5.79) 



Considering Eq. (2.25), we find that Fij now has only two non-vanishing terms, given 
by 

F.i = - (Fi)., + || (Fnhj . (5.80) 

Substituting it into Eqs. (5.73) and (5.74), we obtain 

71 = 1 + 710^1 + ^, g, = 0. (5.81) 


It should be noted that the above solution holds not only in the infrared (IR) regime 
but also in the UV. 

To study the global structure of the spacetime, let us hrst introduce a new radial 
coordinate r* via the relation 


dr 


r = 


1 + ^ 
r 


_ C {Vr + C + y/ry 
~ C 


0, r = 0, 

+ ^/r(r + C) = <( (5.82) 

oo, r = oo. 


In terms of r* the metric takes the form. 


= —dt^ -I- dr*"^ + r'^{r*)d^Vt. 


(5.83) 


Then, one might introduce the two double null coordinates u and v via the relations, 

M = tan“^(t-I-r*), v = — r*), (5.84) 


so that the metric hnally takes the form, 

ds"^ = -5-5-h r^(M, n)d^r2, (5.85) 

cos^ u cos^ V 

where —7r/2 <u,v< 7^/2. The corresponding Penrose diagram is given by Fig. 5.2. 

However, the coordinate transformations (5.84) are not allowed by the foliation 
preserving diffeomorphisms Diff(M, F) of Eq. (2.2). So, in the HL theory the re¬ 
stricted diffeomorphisms do not permit Penrose diagrams. In addition, due to the 
breaking of the general covariance, even if one were allowed to do so, the causal 
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Figure 5.2: The Penrose diagram for = 0, C > 0 and = 0. The double vertical solid 
lines represent the center (r = 0), at which the spacetime is singular. This singularity is 
clearly naked. Note that the restricted diffeomorphisms (2.2) do not allow for the transfor¬ 
mations needed in order to draw Penrose diagrams. Therefore, these diagrams cannot be 
used to study the global structures of spacetimes in the HL theory but are included only 
for comparison. 
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structure of the spacetime cannot be studied in terms of it, as shown explicitly in the 
previous sections for the Newtonian theory. 

Allowed are the coordinate transformations 

t = tanf, r* = tanf*, (5.86) 

where —7r/2 < t < 7r/2 and 0 < f* < 7r/2. Then, the global structure of the spacetime 
is given by Fig. 5.3. 

5.3.1.2. C* < 0. In this case, setting C = —2M < 0, the corresponding metric 
reads, 

X _ ^ 2 ^ 2 ^^ (5.87) 

r J 

This is the solution first found in [73], in which it was argued that the relativistic lapse 
function should he J\f = N — A in the IR. It is not clear how to then relate J\f to N 
and A in other regimes. Instead, in this paper we shall simply take the point of view 
that A and (p are just gravitational gauge helds, and their effects on the spacetime 
itself occur only through the held equations [78]. With the above arguments, we can 
consider the solution valid in any regimes, including the IR and UV. 

Let us hrst note that the metric (5.76) is asymptotically hat and singular at both 
r = 0 and r = 2M. The singularity at r = 0 is a curvature one, as can be seen 
from Eq. (5.78), but the one at r = 2M is more peculiar. In particular, in the region 
r < 2M both t and r are timelike, in contrast to General Relativity where t and r 
exchange their roles across r = 2M. All the above indicate that the nature of the 
singularity at r = 2M now is diherent. In fact, as to be shown explicitly below, the 
region r < 2M actually is not part of the spacetime. 
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Figure 5.3: The global structure of the spacetime in the (i, f*)-plane for N'^ = 0, C > 0 
and Ag = 0. The double vertical solid lines represent the center (r = 0), at which the 
spacetime is singular. The vertical line AB represents the spatial infinity r = oo, while the 
horizontal line i~^A {i~B) is the line where t = oo (t = —oo). The lines t = constant are 
the straight lines parallel to OC, while the ones r = constant are the straight lines parallel 
to The lines BP, BO, BQ, PA,OA and QA represent the radial null geodesics. 
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To see this closely, let us first consider the radial timelike geodesics. It can be 


shown that they are given by, 


t 

T 


Et + to, 

1 

-1 


Min 


(r — M) + \/r{r- — 2M) — y^r{r — 2M) > + ro(5. 


where E is an integration constant, and r denotes the proper time. The constant Tq is 
chosen so that T(ro) = 0 at the initial position of the test particle, r = tq > 2M. The 
sign corresponds to the out-going (in-going) radial geodesics. It is clear 
that, starting at any given hnite radius tq, observers that follow the null geodesics 
will arrive at r = 2M within a hnite proper time. As shown in the last section, 
massless test particles in the HL theory do not follow null geodesics, because of the 
non-relativistic dispersion relations (5.1). In other words, in the HL theory particles 
that follow the null geodesics are not massless and even may not be test particles. 
Setting 

efo) - 1)/, 0, O) , (5.89) 

where / = 1 — 2M/r, we hnd that the spacelike unit vectors. 


efi) = (yE^,-E/7,0,0), 

ef2) = ^(0,0,1,0), 

efs) = (5-90) 

together with form a freely-falling frame, 

efafa (b) = Vab, = 0, (5.91) 


where Da denotes the 4D covariant derivatives, and r]ab is the 4D Minkowski metric 
with a, b, = 0,..., 3. Then, from the geodesic deviations. 


+ fCt’l'' = 0 , 


Dt^ 


(5.92) 


103 



where /Cafe = ~Rcraii'yGla)^^o)^%)^'lb) denotes the tidal forces exerting on the observers, 
we hnd that in the present case /Cafe is given by 

V _ ~ 1)M (^2i;2 I x■2x■3^ nQ^ 

/Cafe ^2 “1“ ' (5.93) 

Clearly, /Cab is finite at r = 2M. All the above considerations indicate that the 
singularity at r = 2M is a coordinate one, and to have a (geodesically) complete 
spacetime, extension beyond this surface is needed. However, unlike that in GR, any 
extension must be restricted to the Diff(M, T) of Eq. (2.2). Otherwise, the resulting 
solutions do not satisfy the field equations. Explicit examples of this kind were given 
in [122]. 

In [73], the isotropic coordinate p was introduced, 

f mV 

r = P f 1 + ^ j , (5.94) 

in terms of which the metric (5.76) takes the form, 

(^dp^ + , (5.95) 

which is non-singular for p > 0. However, this cannot be considered as an extension to 
the region r < 2M, as now the geometrical radius r is still restricted to r G (2M, cx)) 
for p > 0, as shown by Curve (a) in Fig. 5.4. Instead, it connects two asymptotic 
regions, where r = 2M acts as a throat, a situation quite similar to the Einstein- 
Rosen bridge [140]. However, a fundamental difference of the metric (5.95) from the 
corresponding one in General Relativity is that it is not singular for any p G (0, cx)), 
while in General Relativity the metric still has a coordinate singularity at p = M/2 
(or r = 2M) [140]. Therefore, in the HL theory Eq. (5.95) already represents an 
extension of the metric (5.76) beyond the surface r = 2M. Since this extension is 
analytical, it is unique. It is remarkable to note that in this extension the metric has 
the correct signature. 
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Figure 5.4. The function r defined: (a) by Eq. (5.94); and (b) by Eq. (5.101). 


It should be noted that the Einstein-Rosen bridge is not stable in General Rela¬ 
tivity [140]. Therefore, it would be very interesting to know if this is still the case in 
the HL theory. 

To study its global structure, we introduce the coordinate r* by 


1 + ^ dp = Mln ^ 


+P 



{ -oo, p = 0, 
oo, p = oo. 


(5.96) 


Then, in terms of r* the metric can be also cast in the form of Eq. (5.83). Following 
what was done in that case, one can see that the global structure of the spacetime is 
given by Fig. 5.5. 

To compare it with that given in GR, the corresponding Penrose diagram is pre¬ 
sented in Fig. 5.6, although it is forbidden in the HL theory by the foliation-preserving 
diffeomorphisms Diff(M, T) of Eq. (2.2), as mentioned above. 

It is interesting to see which kind of matter helds can give rise to such a spacetime 
in GR. To this purpose, we hrst calculate the corresponding 4-dimensional Einstein 
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C i ^ D 



Figure 5.5; The global structure of the spacetime for = 0, C = —2M < 0 and = 0. 
The vertical line represents the Einstein-Rosen throat (r = = 2M), which is non¬ 

singular and connects the two asymptotically-flat regions I and I'. The horizontal line 
AB (CD) is the line where t = —oo (oo), while the vertical lines CA and DB are the lines 
where r = oo. The lines t = constant are the straight lines parallel to while the ones 
r = constant are the straight lines parallel to The curved dotted lines AD and BC, 

as well as the solid straight lines AD and BC, are the radial null geodesics. 
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Figure 5.6: The Penrose diagram for = 0, C = —2M < 0 and Ag = 0. The straight 
lines represent the future null infinities at which we have r = oo and t = oo, while the 
ones i~i^ represent the past null infinities where r = oo and t = —oo. The vertical line 
represents the Einstein-Rosen throat (r = 2M), which is non-singular and connects the two 
asymptotically-flat regions. 


tensor, 




(5.97) 


which corresponds to an anisotropic fluid, ' 

with = 0, pf^ = M/{AnGr‘^) and = —Mr/{87rG), where Ug = <5* and 
Tg = Clearly, such a fluid does not satisfy any of the energy conditions [129]. 

In particular, when r 3> 1 the tangential pressure becomes unbounded from below, 
while the radial pressure vanishes. Such a fluid is usually considered as non-physical 




in GR. However, in the current setup the spacetime is vacuum, and one cannot elim¬ 


inate it by simply considering the energy conditions. Then, if the configuration is 
stable, one can use it to construct time-machines [141]. 

Inserting Eq. (5.77) into Eq. (5.72), and considering the fact that the range of r 
now is r G {2M,oo), we find that the Hamiltonian constraint is satished, provided 
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that 


A = 0, 20{ge-3gs) -231g3eM^ = 0. 
Then, Eqs. (5.73) and (5.74) have the solution, 

2M g^ 


A — 1 + Aoa/1 


+ 


16 (r — 


r 40C2MV6 L 
8 M^{r + M)r^ - 3M^{5r^ + 7Mr + lOhOM^) 


95 = 98 = 0 . 


(5.98) 


(5.99) 


It is interesting to note that, replacing p by —y we hnd that in terms of y metric 
(5.95) takes the form, 

4 


ds"^ = —dt^ + {dy^ + y^d^Q^ , 


(5.100) 


from which we can see that the geometrical radius now is given by 

-(-T 


(5.101) 


Clearly, the whole region 0 < r < cxd now is mapped to 0 < y < M/2, as shown by 
Curve (b) in Fig. 5.4. Metric (5.100) can be also obtained from metric (5.95) by the 
replacement, M —)■ —M and p ^ y. So, it must correspond to the case C > 0, i.e., 
the one with a negative mass, described in the previous sub-case. 


5.3.2 C = 0, Ag^O 
We have 

for which we hnd that 


u = —- In f 1 — 

2 V 3 ^ 


(5.102) 


r At, ^(^92 + 93) ^2 , 8 ( 9 fi '4 + 35-5 + fi'e) 

Lv - 2 (A - Ag) + -- Ag + - — - Ag, 

^ ^ ^3C"^ (Ag — 3A) -|- 2(‘^ ( 3 g 2 + 93 ) Ag -l- 4 ( 95^4 -|- 3 g 3 + g^) Agj. (5.103) 


F = 

9(4 


To study the solutions further, we consider the cases Ag > 0 and Ag < 0, separately. 
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5.3.2.1. Kg < 0. In this case, defining r-g = y/s/lAj, we find that the corre¬ 
sponding metric takes the form, 


ds^ 



-I- 


(5.104) 


which shows that the metric is not singular except at r = 0. But, it can be shown 
that this is a coordinate singularity. Setting 


r 


* 



(5.105) 


one can cast the metric (5.104) exactly in the form of Eq. (5.83). Then, its global 
structure is that of Fig. 5.3, and the corresponding Penrose diagram is given by Fig. 
5.2, but now the center r = 0 is free of any spacetime singularity. Thus, the range of 
r now is r G [0, cxd). We then find that the Hamiltonian constraint (5.72) is satisfied, 
provided that Cy = 0, i.e.. 


AC\® + 6(3(72 + 93)Crl - 12(9(74 + Sg, + g^) = -3CV|. (5.106) 


Inserting the above into Eqs. (5.73) and (5.74), we obtain the solution. 


A — An\ 1 + 


+ Ai, 


(5.107) 


where Hi is a constant, given by 


Hi = 1 - Ar - 


2 3 — 3g2 — 93 




(5.108) 


5.3.2.2. Ag > 0. In this case, the corresponding metric takes the form. 


d.s'^ = —dt^ + 


dr"^ 


+ r^d^VL. 


1_ il 

'■9 


(5.109) 
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Clearly, the metric has wrong signature in the region r > Vg. In fact, the hypersurface 
r = Tg already represents the geometrical boundary of the spacetime, and any exten¬ 
sion beyond it is not needed. To see this clearly, we first introduce the coordinate r* 
via the relation, 

dr 

f^) 

Then, in terms of r* the corresponding metric can be cast in the form = r'^ds'^, 
where 

ds"^ =—dt^ + dx"^ + sin^ xd'^il, (5.111) 

with t = t/vg, X = r*/vg. But, this is exactly the homogeneous and isotropic Einstein 
static universe, which is geodesically complete for—cxD<t<cxD, 0<a;<7 r, 0< 
9 < 71 and 0 < 0 < 27i, with an i? x topology [129]. Then, it is easy to see that its 
global structure is given by Fig. 5.3, but now the vertical line is free of spacetime 
singularity, and the line AB is the one where r = Vg (or x = ti). The corresponding 
Penrose diagram is given by Fig. 5.7. 

Therefore, in this case the range of r is r G [0, r^j. Then, the Hamiltonian 
constraint (5.72) requires. 



= Tg arcsm 


(5.110) 


+ 6 ( 35^2 + + 12( 95^4 -|- ?>g^ -|- — 3C^’"g- (5.112) 

Hence, Eqs. (5.73) and (5.74) have the solution, 

A = Ao^l- +H 2 , (5.113) 

where A 2 is another integration constant, given by 

^ 9 

It should be noted that in General Relativity the Einstein static universe is ob¬ 
tained by the exact balance between the gravitational attraction of matter (p^ = 
Pc Pm = 0) and the cosmic repulsion (A = Ac), where Ac = AnGpc- As a result, the 

no 



• + 



Figure 5.7: The Penrose diagram for = 0, (7 = 0 and > 0, which is the Einstein 
static universe. The curves i~i^ and i~^i^ are, respectively, the lines where t = —oo, x = tt, 
and t = + 00 , x = n. 
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configuration is not stable against small perturbations [142]. However, in the present 
case since the spacetime is vacuum, Eq. (5.112) suggests that the balance is made 
by the attraction of the high-order curvature derivatives and the cosmic repulsion, 
produced by both A and A^. Then, it would be very interesting to know whether it 
is stable or not in the current setup. 


5.3.3 C ^0, Ag^O 

When Ag, C 7 ^ 0, we hnd that 

9C^ + SA^r® 


= 


6 r® 


1 

D* DJ _ 

n^nkUi - 


27C'^ + lOSAgCV^ + 32A^r® 




(5.115) 


from which one can see that the spacetime is singular at r = 0. Moreover, we hnd 
from (2.25) that 

1 


36r8C^F(r) 


-27^^(22(75 + 25(76 - 20(78) - 81C'V(8(75 + 9g^ - 7gs) 


-9C\^ [Ag{-2Qg^ - 30(76 + 25 ( 78 ) + 

-|-12Cr® —3C^ -|- Ag^^( 12(72 -|- hg^) + A^(36(74 -|- 14(75 + 6(75 
-|-4r® [—3C‘^(3A — Ag) + 2(‘^AR3g2 + ( 73 ) + AARdg^ + 3g^ -|- g^) I /■(,5.116) 


where the third-order polynomial F{r) is dehned by 

A. 


F{r) = C + r - -r^ 


i.e. 




F(r) 


The function Cy is given by 


Cy — 


a 


-f- /Sr -f -|- 5r‘^ 


36r^C 


(5.117) 
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where 


a = 27C'3 c/ 6 + SlOC'^^s, 

(3 = 810 ^ 2 ( 78 , 

7 = lOSC'^^^sAg + 108C'2(76A5 - 270C'2^8A5 + 54C'2^/3C^ 

<5 = 14%2C'A2 + 288(74A^ + 96(/5A^ + 32(/6A^ + 48(73C'A2 + 72C'A-72CX- 

All the quantities in (5.115) are hnite for any r 7 ^ 0. On the other hand, from 
Eq. (5.71) one can see that the metric coefficient Qrr could become singular at some 
points. To study the nature of these singularities, we distinguish the four cases, 
O > 0 , Ag > 0 ; O > 0 , Ag < 0 ; O < 0 , Ag > 0 ; and O < 0 , A^ < 0 . 


5.3.3.1. 0 > 0, Ag > 0. In this case, the polynomial F{r) has exactly one real 
positive root at, say, r = rg{C,Kg) > 0, as shown in Fig. 5.8. We hnd that 


e 


2v 


r 

D{r){rg — r) ’ 


(5.118) 


where D{r) = Ag{r‘^ + r^r + d)/3, d = — 3/Ag, and D{r) > 0 for all r > 0. 

Introducing the coordinate x via the relation 


/ dv 

v?=^ = (5.119) 

or, inversely, r = Vg — x'^, the corresponding metric in terms of x takes the form 

ds'^ = —dt^ + ^ d^x + (r„ — x^Yd^Vt, (5.120) 

D[x) ^ ' 

where D{x) = Ag(a;"^ —+ —3/Ag)/3 > 0 for |a;| < Clearly, the coordinate 

singularity at r = (or a; = 0 ) now is removed, and the metric is well dehned for 
I a; I < At the points, x = ± 1/^7 (or r = 0 ), the spacetime is singular, as shown 
by Eq. (5.115). Thus, in the present case the spacetime is restricted to the region 
|a;| < —00 < t < cxD in the (t, a;)-plane, with the two spacetime singularities 
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located at a; = as its boundaries. The global structure of the spacetime and 

the corresponding Penrose diagram are shown in Fig. 5.9. 

The change of variables (5.119) can be understood by considering the one-form 

dr = 

Even though the denominator of the right-hand side vanishes at r = r^, we can turn 
e‘'dr into a nonsingular one-form by introducing a Riemann surface. Indeed, if we 
promote r to a complex variable and dehne the genus 1 Riemann surface S as the two- 
sheeted cover of the complex r-plane obtained by introducing two branch cuts along 
the intervals [0,rg] and [ri,r 2 ], where ri and r 2 are the two (possibly complex) zeros 
of D{r), e'^dr is a holomorphic one-form on S. Letting (0,rg]i and (0,rg]2 denote 
the covers of the interval (0,rg] in the first and second sheets of S, respectively, 
the spacetime consists of points (r, 6*,0, t) with r G (0,rg]i U (0,rg]2. The variable 
X = — r introduced in (5.119) is analytic near the branch point at r = and 

r G (0,rg]i U (0,rg]2 corresponds to a; G (— ^/^)■ We can fix the dehnition of x 
by choosing the branch of the square root so that, say, a; > 0 for r G (0,rg]i. Thus, 
in terms of the variable a;, the spacetime manifold can be covered by a single global 
chart (no double cover is necessary) and the metric ds‘^, which involves the square of 
the differential e‘'dr, is manifestly nonsingular at r = Vg. In particular, the metric of 
the extended spacetime is analytic, which ensures that the extension is unique. 

The Hamiltonian constraint is 

Cve^r'^dr = 0. (5.121) 

Indeed, the Hamiltonian constraint (2.15) should be interpreted as 

j Cv^o\g = 0, (5.122) 

where Vol^ is the volume form induced by the metric Qij and the integration extends 
over a spatial slice of the spacetime. Using the variables (r, 9, 0), we have 

Volg = sin OdrdOdcj), 
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Figure 5.8: The function F{r) = re for W = 0, C > 0 and Ag > 0, where 



Figure 5.9: (a) The spacetime in the {t, x)-plane, where xq = (b) The Penrose 

diagram for N'^ = 0, C > 0, > 0. The curves are the lines where r = 0, at which 

the spacetime is singular. The straight line i~i~^ represents the surface r = rg. 
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and the integration extends over 9 G [0,7r], 0 G [0,27r], and r G [0,rg]i U [0,rg]2. By 
symmetry, the contributions from the sets where r G [0, and r G [0, rg ]2 are equal. 
Since each contribution is proportional to the left-hand side of (5.121), the constraint 
reduces to (5.121). 

In view of (5.117), the constraint (5.121) becomes 

a + l3r + -|- 5r^ 


-\/rdr = 0. 


(5.123) 


Jo 36r'^('^^/F{r) 

Denoting the integrand in (5.123) by /(r), we see that |/(r)| is bounded by a constant 
times 1 / \/i^g — f as r ^ Vg. Thus, the integral converges near Vg. On the other hand, 
as r —)■ 0, 


/(r) = 


a 


+ 0 




so that (5.123) can only be satished if a = 0. Using similar arguments, we infer that 
the coefficients /9,7, S must also vanish, i.e., 

a = f3 = 'y = 5 = 0. 

Solving these equations, we conclude that the Hamiltonian constraint is satished if 
and only if the g/s satisfy the following four conditions: 

C^(A, - A) - 2g2eAl 


9 a = 


9b = 


4A3 


(5.124) 


2A/ 


96 = 0, 98 = 0. 


Using the conditions (5.124) in the expression (5.116) for F^r, we hnd that Eqs. (5.73) 
and (5.74) have the solution 

'/W) r FrriF){FfAdF 


A{r) = 


2^- 


'ro 


\fW) 


(5.125) 


where 


F = 

J. rfy. 


- 2970^gs - 32Aa^rgo + 126CVAggs 


36r^(‘^AgF{r 

+ 12Cr^ \2Al{3g2 + go) + C'(9A - 6A,) 


+8r''A„ 


2A;(3(72 + (73) + C'(9A-6A,) 


(5.126) 
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and ro G (0,rg) is a constant. The integrand in (5.125) is smooth for 0 < r < r^. 
Thus, A{r) is a smooth function of r G (0,rg). Unless = 0, the integral diverges 
as r —)■ 0, so that A{r) has a singularity at r = 0. As r —)■ r^, the integrand is 
bounded by const x {vg — This implies that A{r) is bounded as r —)■ Vg. In 

fact, viewed as a function on the Riemann surface S, A{r) is analytic near r = Vg. 
This follows since the integrand in (5.125) is a meromorphic one-form with a pole of 
at most second order at r = r^. Thus, the integral has a pole of at most order one at 
Tg, which is cancelled by the simple zero of the prefactor \jF{r) = D{r){r-g — r). 
In conclusion, the gauge held A given by (5.125) is a smooth function everywhere on 
the extended spacetime away from the singularity at r = 0. 

5.3.3.2. C > 0, Ag < 0. In this case, F(r) > 0 for r > 0 and the metric 
coefficient Qrr is positive and non-singular except at the point r = 0, at which a 
naked spacetime singularity appears. The corresponding Penrose diagram is given by 
Fig. 5.2 with r G (0,cxd). The Hamiltonian constraint (5.72) requires that 

Cve’^r^dr = 0. (5.127) 

As in the previous subsection, this constraint is equivalent to the conditions given in 
(5.124). 

The function A{r) is again given by the formulas (5.125)-(5.126) and is a smooth 
function of r G (0, cxd). As r — )■ cx), the absolute value of the integrand is bounded by 
constant x r“^. Thus, choosing ro = cxd in (5.125), we hnd that A{r) is bounded as 
r —)■ CX). Unless = 0, the integral diverges as r —)■ 0, so that A{r) has a singularity 
at r = 0. 



5.3.3.3. U < 0, Ag > 0. In this case, if Ag > 4/(9U^), = r/F{r) is strictly 

negative for all r > 0, so that, in addition to t, the coordinate r is also timelike. The 

physics of such a spacetime is unclear, if there is any. Therefore, in the following we 
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consider only the case 


0 < A, < ^. (5.128) 

Then, we hnd that F{r) is positive only for 0 < r_ < r < r+, where r±{Kg,C) are 
the two positive roots of F{r) = 0, as shown in Fig. 5.10. We write as 




(5.129) 


(r + ro)(r — r_)(r+ — r) ’ 

where ro(Ag,C') > 0. To extend the solution beyond r = r±, we shall hrst consider 
the extension beyond r = r_. Such an extension can be obtained via 


= yjr-r_, (5.130) 

or inversely, r = + r_. Since r < r+, we hnd that —xq < x < xq with xq = 

— r_. It can be seen that the coordinate singularity at r = r_ disappears, and 
the extended region is given by |a;| < Xq, as shown by Fig. 5.11 (a). 

To extend the solution beyond r+, we introduce x via the relation 



r = r+ - (x =F XqY, 


(5.131) 


where the ” sign applies when x > Xq and the “+” sign applies when x < —Xq. 
Fig. 5.11 (b) shows the graph of r as a function of x. From Fig. 5.11 we can see 
that the extension along both the positive and the negative directions of x need to 
continue in order to have a maximal spacetime. This can be done by repeating the 
above process inhnitely many times, so hnally the whole (t, a;)-plane is covered by an 
inhnite number of hnite strips, in each of which we have r_ < r < r_|_. The global 
structure is that of Fig. 5.12 and the corresponding Penrose diagram is given by Fig. 
5.13. Thus, in this case we have r G [r_,r_|_]. 

The Hamiltonian constraint (5.72) requires that 



a + /3r + 

36F('^^/F(r) 


^/rdr = 0 . 


(5.132) 
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Figure 5.10: The function F{r) = re~‘^'' for C < 0 and Kg > 0, where 

F{r) = 0 has two positive roots r± only for Kg < 4/(9C'^). When Kg > 47(9(7^), F{r) is 

always non-positive for any r > 0. 



(a) (b) 

Figure 5.11: (a) The function r vs x given by Eq. (5.130), where xq = \/r+ — r_. (b) The 
function r vs x given by Eq. (5.131), where xi = xq + \/To- 
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E A i^B F 



Figure 5.12: The global structure of the spacetime for C* < 0, > 0 and Kg < 4/(9C'^). 

The vertical line is the one where r = r-, and the ones AC and BD represent the 
lines where r = r_|_, while on the lines EG and FH we have r = r-. The spacetime repeats 
itself infinitely many times in both directions of the x-axis. 



Figure 5.13. The Penrose diagram for C < 0, Kg > 0 and Kg < 47(9(7^). 
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Geometrically, this condition can be understood by introducing a Riemann surface 
S as a double cover of the complex r-plane with two branch cuts along [r_,r+] and 
[—ro,0]. The integrand in (5.132) is a one-form a; on S which is holomorphic in 
a neighborhood of the closed curve ai = [r_,r+]i U [r+,r_] 2 . Topologically, the 
elliptic curve S is a torus, ai is a nontrivial cycle, and the condition (5.132) states 
that the integral of u along the cycle oi vanishes. This imposes a constraint on the 
coefficients which translates into a condition on the g/s involving elliptic 

integrals. Assuming this condition to hold, the function A{r) is given by (5.125) with 
tq G (r_,r+) and Frr as in (5.116). 


5.3.3.4- C < 0, Ag < 0. In this case, the function F(r) = is positive only 
for r > Tg, as shown in Fig. 5.14. Thus, can be written in the form, 

r 




(5.133) 


D{r){r — Tg) ’ 

where D{r) > 0 for r > 0. The extension can be carried out by introducing a new 
coordinate x via the relation. 


r = r-g. (5.134) 

In terms of x the coordinate singularity at r = disappears, and the extended 

spacetime is given by —oo < t,x < oo m the (f, a;)-plane. Its global structure is given 

by Fig. 5.5, while the corresponding Penrose diagram is given by Fig. 5.6. Thus, in 

this case the range of r is r G [r^, cx)). 

The Hamiltonian constraint (5.72) requires that 

a + (3r + ^ 

/ - , —yrdr = U. 

Jr, 

The behavior of the integrand as r —)■ cxd implies that 


q; = /9 = 7 = (5 = 0, 

so that the constraint reduces to (5.124) and the function A{r) is given by (5.125)- 

(5.126), which is not singular everywhere in the extended spacetime. 
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root of F{r) = 0. 


5.4 Vacuum Solutions with N'' 7 ^ 0 
When 7 ^ 0, the vacuum solutions are given by [78], 

ds^ = —dt^ + [dr + e'^~'^dtY + r'^d^Vt, 


with 


j. , , 2m 1,2 ^ / \ r A / 1 ! \ 

fjt, = 2 I -^ ~ 2A{r) H— / A{r )dr ), 


u = (p = Ag = 0, 


(5.135) 


(5.136) 


where the gauge held A must satisfy the Hamiltonian constraint, 


rA'{r)dr = 0. (5.137) 

Otherwise, it is free. However, as shown in [78], the solar system tests seem uniquely 
to choose the Schwarzschild solution A = 0. Therefore, in the following we shall 
consider only this case. 



/i 



z/ 


(p = Ag = A = 0. 
122 


(5.138) 



It should be noted that if (iV, z/, N^) is a solution of the vacuum equations, so is 
(TV, z/, — TV’’). The latter can be easily obtained by the replacement t —)■ —t. With 
such changes, we have Kij —)■ —Kij (in the static case). Clearly, these do not affect 
the singularity behavior. We then obtain [122, 103] ^ , 


Rij — 

K = 


0 , 

Cl 


]j r^ ( 6 m + Ar») 0 “ +''C . 




27m? + 6 mAr^ + A^r® 


v.."? (a 


(5.139) 


where ei(= ±1) originates from the expression TV’’ = eic^, obtained by the replace¬ 
ment t —t, as mentioned above. To further study the above solutions, let us 
consider the cases (1) m = 0 ,A 7 ^ 0; ( 2 ) m 7 ^ 0 ,A = 0; and (3) m 7 ^ 0,A 7 ^ 0 
separately. We shall assume that m > 0, while A can take any values. 


5.4.1 m = 0 , A^O 

In this case, only A > 0 is allowed [78], as can be seen from Eq. (5.136). That 
implies that the anti-de Sitter spacetime cannot be written in the static form of 
Eq. (5.5) with the projectability condition. Then, we have TV^ = / = 1, TV’’ = eir/t', 
or 

ds'^ = —dt^ + (^dr + + r‘^d?Vt, (5.140) 

where 4 = a/3/|A|. Without loss of generality, we shall consider only the case ei = —1, 
as the case ei = 1 can be simply obtained from the one ei = — 1 by inverting the time 
coordinate. In terms of TV, TV*, gij or their inverses, TVj, the metric is non-singular, 
except for the trivial r = 0 and 6* = 0,7r. In addition, from Eq. (5.139) we also hnd 
that 

K = -7^, = A, (m = 0). (5.141) 

^ There is a typo in the expression of K given by Eq. (3.2) in [122]. Although it propagates to 
other places, this does not affect our main conclusions, as K and KijK^^ have similar singularity 
behavior. 
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On the other hand, in terms of the 4-dimensional metric, and , it is not singular 
either, as one can see from the expressions. 





(5.142) 


although the nature of the radial coordinate does change, 

timelike, r > i, 

null, r = i, (5.143) 

spacelike, r < i. 

To study the solution further in the HL theory, we consider two different regimes, 
E and E 3> M*, where M* = min. {Ma, Mb, ...} and M„’s are the energy 

scales appearing in the dispersion relation (5.59). 



5.4-1-1. E ^ M*. When the energy E of the test particle is much less than M*, 
from Eq. (5.59) one can see that E{() ~ (. This corresponds to the relativistic case 
(n = 1). Then, for the ingoing test particles (e = —1), we have 

H = N^+N^ = ^. (5.144) 

Thus, the hypersurface r = i is indeed a horizon. In fact, it represents a cosmological 
horizon, as first found in General Relativity [143]. 

However, because of the restricted diffeomorphisms (2.2), it is very interesting to 
see the global structure of the de Sitter spacetime in the HL theory. To this purpose, 
let us consider the coordinate transformations, 

t’ = r' = (5.145) 
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in terms of which the corresponding metric takes the form, 


ds"^ = -dt^ + [dr'‘^ + . 

(^—dt'"^ + dr'"^ + r'‘^d^Q) . (5.146) 

From Eq. (5.145) we can see that the whole (t,r)-plane, —oo < t < oo, r > 0, 
is mapped to the region t', r' > 0. However, the metric now becomes singular at 
F = 0, cxD (or t = ± cxd ). To see the nature of these singularities, one may recall the 
5-dimensional embedding of the de Sitter spacetime in General Relativity [129], from 
which we hnd that in terms of the 5-dimensional coordinates v and w, t' is given by 
t' = /{v+w). Therefore, F > 0 corresponds to v+w > 0. Thus, the region F, r' > 0 

only represents the half hyperboloid n -f tc > 0, as shown by Fig. 16 (ii) in [129]. 
In particular, t' = 0 represents the boundary of the spacelike inhnity, so extension 
beyond this surface may not be needed. Although the extension given in [129] in 
terms of the static Einstein universe coordinates (t, y, d, 0) is forbidden here by the 
restricted diffeomorphisms (2.2), as that extension requires, 

cos(x) + sinh 

the extension across t' = oo (or n -f- tc = O’*") seems necessary. 

Another way to see the need of an extension beyond = 0 is that the metric 
(5.146) is well-dehned for t' < 0. So, one may simply take —oo < t' < oo. But, 
this cannot be considered as an extension, as the metric (5.146) is singular at t' = 0, 
and the two regions t' > 0 and t' < 0 are not smoothly connected in the t',r'- 
coordinates. In this sense, a proper extension is still needed. However, due to the 
restricted diffeomorphisms (2.2), it is not clear if such extensions exist or not. Fig. 
5.15 shows the global structure of the region t' > 0, which is quite different from its 
corresponding Penrose diagram [143]. 



t = f'ln 


cosh ( - 
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A B 



Figure 5.15: The global structure of the de Sitter solution = / = 1, N'' = —\frjl in the 
HL theory with the restricted diffeomorphisms (2.2) for the region t' > 0. The horizontal 
line AB corresponds to = oo (or t = —oo), while the vertical line BD to r' = oo (or 
r = oo). 


5.4-1.2. E ^ M*. When the energy E of the test particle is greater than M*, 
from Eq. (5.59) one can see that high order momentum terms become important, and 
E{() ~ {n >2). For the sake of simplicity, we consider the case with n = 2 only. 
Then, from Eqs. (5.25) and (5.30) we hnd that 


X 

H 


2iE 

Vr^ + APE^' 

r AiE 

i 1 / 7-2 _|_ 4PE + r 


(5.147) 


Thus, H{r, E) = 0 has only one real root, 

at which we hnd that 


th = 


H{rH,E) = 


12iE 


> 0 . 


(5.148) 


(5.149) 


Ai‘^E + 3rjj 

Eqs. (5.39) and (5.40) then tell us that the surface r = th is a horizon for a test 


particle with energy E. It should be noted that, in contrast to the Schwarzschild case 
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[cf. Eq. (5.44)], th now is proportional to E, that is, the higher the energy of the test 
particle, the lager the radius of the horizon. To understand this, let us consider the 
acceleration of a test particle with its four-velocity u\ = —5*, located on a surface r. 
Then, we find that 

{ —^5]], Schwarzschild, 

(5.150) 

deSitter. 

That is, for the Schwarzschild solution, the test particle feels an attractive force, while 
for the de Sitter solution, it feels a repulsive one. Because of this difference, in the de 
Sitter spacetime th is proportional to E, in contrast to the Schwarzschild one, where 
it is inversely proportional to E, as shown explicitly in Eq. (5.44). 


5.4-.2 m > 0, A = 0 

When A = 0 and m > 0, it is the Schwarzschild solution. In particular, in the IR, 
the surface r = 2 m represents a horizon, while for high energy particles, the radius of 
the horizon is energy-dependent, as explicitly given by Eq. (5.44) for n = 2. So, we 
shall not repeat these studies, but simply note that now the solution takes the form, 

ds'^ = —dt^ + — \j -^dt 

which is singular only at r = 0, as can be seen from Eq. (5.139). So, it already 
represents a maximal spacetime in the HL theory. 

It is interesting to note that the above metric covers only half of the maximally 
extended spacetime given in GR. This can be seen easily by introducing the coordinate 
r [103], 


j (5.151) 


r 


t — 




mr 


dr 


r — 2m 
t — 2\/2mr — 2m In 


r — 2 m 


(yT -h y/^y 


(5.152) 
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in terms of which, the solution takes the standard Schwarzschild form, = —/(r)dr^+ 
f~^{r)dr‘^ + r'^dVL^ with /(r) = 1 — 2m/r. Of course, the above transformations are 
forbidden by Eq. (2.2). 

5 . 4.3 m>0, A^O 

In this case, it is convenient to further distinguish the two subcases A > 0 and 
A < 0. 


5. 4 .3.1. A > 0. In this case, the metric takes the form, 

2 


ds^ = —dr + \ dr — 


r r 


(5.153) 


When E M*, as in the last case the dispersion relation becomes relativistic, 
and ~ which we have n = 1. Then, we find that 


H{r) = l + + ^ 

F(r) 


(1 + / 


2m r 


(5.154) 


but now F{r) = — (r^ — ^‘^r + 2m(?). Clearly, F{r) has one maximum and one mini¬ 
mum, respectively, at r = where = ^/\/3 and F{rm) = —2i‘^{m — 1/(3\/A), 
as shown in Fig. 5.10. Thus, when > 1/(9A^), H{r) = 0 has no real positive root, 
and a horizon does not exist even in the IR. Therefore, the singularity at r = 0 is 
naked. When < 1/(9A^), H{r) = 0 has two real and positive roots, r±, (r+ > r_), 
where r = r+ is often referred to as the cosmological horizon and r = r_ the black 
hole event horizon [143]. When = 1/(9A^), the two horizons coincide. In GR, 
the corresponding Penrose diagrams were given in [143]. However, as argued above, 
in the HL theory these diagrams are not allowed, as they are obtained by coordi¬ 
nate transformations that violate the restricted diffeomorphisms (2.2). Nevertheless, 
since the metric is not singular in the current form, it already represents a maximal 
spacetime. 
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When E 3> M*, the high momentum terms dominate, and for n = 2, we hnd that 


where A(r) > 0 for any r G (0, cxd), and F{r) = — 4Ei'^r/3 + 2m£‘^. It can be 

shown that when > S£E^^'^/27, H{r) = 0 has no real and positive roots. Thus, in 
this case there are no horizons, and the singularity at r = 0 must be naked. When 
m? < 8£E^/‘^/27, H{r) = 0 has two real and positive roots, say, ri ^2 {^2 > ri), but 
now ri ^2 = 'Yi^2{E,Yn, i). Thus, in this case there also exists two horizons, but each of 
them depends on E. When m? = 8f'i?^/^/27, we have ri = r 2 , and the two horizons 
coincide. 


5.4-3.2. A < 0. In this case, the metric takes the form, 

ds'^ = —dt^ + ^ + r^dVt^, (5.156) 

where £ = a/3/|A|. Then, from Eq. (5.139), it can be seen that the spacetime 
is singular at = {2m£^Yf^ [122]. This is different from GR, in which the only 
singularity of the anti-de Sitter Schwarzschild solution is at r = 0. 

When E M*, as in the last case the dispersion relation becomes relativistic. 


Then, we find that 


H{r) 


1 + = 1 - 
E{r) 



j.£2 



£2 


(5.157) 


but now with E{r) = + £‘^r — 2m£‘^, which is a monotonically increasing function, 

as shown by Fig. 5.14. Thus, H{r) = 0 has one and only one real and positive root 
= fHi'rn,£). But, is always less than r^, i.e., vh < r^. Thus, the singularity at 
r = r* is a naked singularity. 
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When E 3> M*, let us consider only the case n = 2. Then, we hnd that 


A'(r) = 


2E 




TT, \ F(r) 

H{r) = \ -— — 2X=———, 

' r A(r) 


(5.158) 


where A(r) > 0 for any r G (0, cxd), and F{r) = + 4Ei‘^r/3 — 2mi‘^. It can be 

shown that this F{r) is also a monotonically increasing function, as shown by Fig. 
5.14, and F{r) = 0 has only one real and positive root, th = rnim, E,i). Again, 
since H{rs) = 1 and H{rH) = 0, we find that th is also always less than r^, although 
now r'H depends on E. Thus, the singularity at r = r* is a naked singularity. 


5.5 Slowly Rotating Vacuum Solutions 

Slowly rotating vacuum solutions in other versions of the HL theory have been 
studied by several authors [144]. The goal of this section is to derive slowly rotating 
black hole solutions in the HMT setup. We will seek a solution of the form 

+ r^(dd^ + sin^ dd0^) (5.159) 

[dr + — au{r) sin^ dd0)]^ , 

where the functions and u{r) are independent of (t,9,(j)). By requiring 

that the metric satisfy the equations to first order in the small rotation parameter a, 
we will be able to determine z/, /i, and u. 

The ansatz (5.159) is motivated by the fact that it agrees with the Kerr solution 
to first order in a. Indeed, the Kerr line element expressed in Doran coordinates [145] 
is given by 




Kerr 


—dt^ + (r^ + cos^ 6)d6‘^ 
+ (r^ + a^) sin^ dd0^ + 


(5.160) 


X 


2 , I r^ + g^cos^d 

+ 0 ? 

\/2mr^r2'+'a^ . 2 


dr + 


+ a? cos^ 9 


-{dt — a sin^ 9d(j)) 
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where m and a are parameters. As a —)■ 0, this metric coincides with (5.159) to first 
order in the rotation parameter a, provided that 


z/(r) = 0, /i(r) = log Y-, a;(r) = 1. 

In particular, when a = 0, it reduces to the Schwarzschild metric in PainlevA 
Gullstrand form. 

Note that the form (5.159) of the line element is compatible with the projectability 
condition N = N{t)] its ADM coefficients are 

N = l, AT* = 


Working in the gauge (p = 0, the momentum constraint (2.17) for the metric (5.159) 
reduces to 


r 


+ 0(a2) = 0, 


2r4 

-2a; (l 


■ (ca" + u' (V 

— r^/i" + - 


2r^/i'^ + rz/') 


+ O(a^) — 0, 


(5.161) 


while the equation (2.20) obtained from variation with respect to A yields 

(1 - r^Ag) e^‘' + 2ru' - 1 + 0{a^) = 0. (5.162) 

The first equation in (5.161) implies that u is constant, and then (5.162) shows that 


z/ = 0, Ag = 0. 


This yields 


R.. = 0(a2), (5.163) 

jCk = + 2rp') + O(a^), 

Cv = 2A+ 0{a^). 
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The (rr)-component of the dynamical equations (2.22) gives 

2rA q — r^A + 2re^^fj,' + ^ 2A'^ 

lyt 

+ 0{a^) = 0, 

where we have assumed that A{r) has the form 

A{r) = Ao(r) + Ai(r)a + 0{a^). 


(5.164) 


(5.165) 


The terms of 0{1) in (5.164) imply that 

/i(r) = ^ In - 2Ao(r) + ^ ^ ^o(s) ds ) , 


(5.166) 


where ro > 0 is a constant, while the terms of 0{a) imply that Ai is a constant. With 
these choices, all the components of the dynamical equations as well as the equations 
obtained from variation with respect to A and (p are satished to hrst order in a, and 
the Hamiltonian constraint (2.15) becomes 

poo 

/ rAQ{r)dr + O(a^) = 0. 

Jo 

Finally, the second equation in the momentum constraint (5.161) is satished to 0{a) 
provided that 

' di 


uj(r) = e 


+ d2r^ 


(5.167) 


di + d2r^ 


2m+ 2 Ao{s)ds - 2rAo + ’ 

where di and d 2 are the integration constants. 

In summary, the ansatz (5.159) gives a solution to hrst order in a provided that 
/i(r) is given by (5.166), u{r) is given by (5.167), and 


u — 0, A(r) — Ao('r) + aAi + C7(a^), 


(5.168) 


where ro > 0 ,m,A,Ai,di,d 2 are arbitrary constants and Ao(r) can be freely chosen 
as long as 

rA'Q{r)dr = 0 . 
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We recover the slowly rotating version of the Kerr solution by taking Aq = 0, A = 0, 
and d 2 = 0. Setting a = 0, on the other hand, we recover the static solutions obtained 
in [78], 

Let us point out that the standard Einstein equations also allow for a nonzero 
value of d 2 in the slowly rotating limit. Indeed, substituting the ansatz (5.159) with 
z/ = 0 into the vacuum Einstein equations 

Rai3-^gaf^R = 0, a,/3 = 0,1,2,3, 


we hnd that they are satished to order 0{a) if and only if 



where m > 0 is a constant, and u{r) is given by (5.167) with arbitrary constants di 


and d 2 - 


5.6 Conclusions 

In this chapter, we have systematically studied black holes in the HL theory, using 
the kinematic method of test particles provided by KK in [135], in which a horizon 
is dehned as the surface at which massless test particles are inhnitely redshifted. 
Because of the nonrelativistic dispersion relations (5.1), we have shown explicitly the 
difference between black holes dehned in General Relativity and the ones dehned 
here. In particular, the radius of the horizon usually depends on the energy of the 
test particles. 

When applying this dehnition to the spherically symmetric and static vacuum 
solutions found recently in [73, 79, 78], we have found that for test particles with 
sufficiently high energy, the radius of the horizon can be made arbitrarily small, 
although the singularities at the center can be seen in principle only by test particles 
with inhnitely high energy. 

We paid particular attention to the global structures of the static solutions. Be¬ 
cause of the restricted diheomorphisms (2.4), they are dramatically diherent from 
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the corresponding ones given in GR, even the solutions are the same. In particular, 
the restricted diffeomorphisms (2.4) do not allow us to draw Penrose diagrams, al¬ 
though one can create something similar to them, for example, see Figs. 5.3, 5.5, 
5.9, 5.12, 5.15. But, it must be noted that, since the speed of the test particles in 
the HL theory can be inhnitely large, the causality in this theory is also dramatically 
different from that of General Relativity [cf. Fig.5.1]. In particular, the light-cone 
structure in General Relativity does not apply to the HL theory. Among the static 
solutions, a very interesting case is the one given by Fig. 5.5, which corresponds to an 
Einstein-Rosen bridge. In GR, this solution is made of an exotic fluid as one can see 
from Eq. (5.97), which is clearly unphysical, and most likely unstable, too. However, 
in the HL theory, the solution is a vacuum one, and it would be very interesting to 
see if this conhguration is stable or not in the HMT setup. 

Finally, we have studied the slowly rotating solutions in the HMT setup [73], and 
found explicitly all such solutions, which are characterized by an arbitrary function 
Ao(^")- When Aq = 0 they reduce to the slowly rotating Kerr solution obtained in 
GR. When the rotation is switched off, they reduce to the static solutions obtained 
in [78]. 
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CHAPTER SIX 


Summary and Outlook 
6.1 Summary 

Isaac Newton was the first person to give a mathematical theory of gravity, which 
he hrst published in his Philosophice Naturalis Principia Mathematica on July 5th, 
1687 [146]. It successfully predicted the return of Halley’s comet, and later many 
planets were discovered using this law of gravitation. Although there were a few 
inconsistencies found in the theory toward the end of nineteenth century, no one 
dared to suggest that a better theory of gravitation was required to account for those 
discrepancies. After the success of the Special Theory of Relativity, Albert Einstein 
saw that Newtonian gravity contradicts relativity. After ten years of work using the 
best mathematical tools of the day, Einstein formulated General Relativity in its hnal 
form a century ago on November 25th, 1915. At that time and for many years later, 
Einstein’s theory was better than Newton’s only to account for three very tiny effects. 

Einstein’s Classical Field Theory is the best description of gravity we have today 
and it describes the universe at large. However, the structure of bodies themselves are 
held together by other fundamental interactions in nature, which are best described 
by Quantum Field Theory. The hnal stages of gravitational collapse and earliest 
moments of the birth of the universe require us to have the same description of all 
the fundamental interactions. This is a daunting task and our current best hope 
is String Theory. The heart of the problem is to formulate gravity as a quantum 
held theory and this problem is often simply called qnantnm gravity. Althongh at 
snch high energies one cannot separate pnre gravity contribntions from that of other 
interactions, it is still a big step forward just to quantize pnre gravity. Many such 
ehorts in this direction are Loop Qnantnm Gravity, Causal Dynamical Triangulations, 
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Asymptotic Safety and the latest Hofava-Lifshitz gravity which is investigated in this 
dissertation. 

Hofava-Lifshitz Theory of qnantnm gravity is motivated by the Lifshitz theory 
in solid state physics. One of the essential ingredients of the theory is the incln- 
sion of higher-dimensional spatial derivative operators which dominate in the nltra- 
violet, making the theory power-connting renormalizable. The exclusion of higher¬ 
dimensional time derivative operators, on the other hand, guarantees that the theory 
is unitary. The problem of non-unitarity has plagued the quantization of gravity for 
a long time. However, this asymmetrical treatment of the space and time variables 
inevitably leads to the breaking of Lorentz and Diffeomorphism symmetries. 

We studied spherically symmetric static spacetimes generally hlled with anisotropic 
fluid. We found the vacuum solutions and used them to understand the solar system 
tests. We found that the theory is consistent with solar system tests, by properly 
choosing the parameter involved in the theory. We considered the gravitational col¬ 
lapse of spherical fluid and derived the junction conditions across the collapsing star 
under the assumption that they be mathematically meaningful in terms of distri¬ 
bution theory. We systematically studied black holes in the Hofava-Lifshitz Theory 
by following the kinematic approach, in which a horizon is dehned as the surface at 
which massless test particles are inhnitely redshifted. Because of the nonrelativistic 
dispersion relations, the speed of particles is unlimited, and test particles do not fol¬ 
low geodesics. As a result, there are signihcant differences in causal structures and 
black holes between General Relativity and HL theory. In particular, the horizon 
radii generically depend on the energies of test particles. By applying them to the 
spherical static vacuum solutions, we hnd that, for test particles with sufficiently 
high energy, the radius of the horizon can be made as small as desired, although the 
singularities can be seen in principle only by observers with inhnitely high energy. 
In these studies, we pay particular attention to the global structure of the solutions. 
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and find that, because of the foliation-preserving diffeomorphism symmetry, they are 
quite different from the corresponding ones given in GR, even though the solutions 
are the same. 


6.2 Outlook 

So far the studies show that Hofava-Lifshitz Theory is not ruled out by any ex¬ 
periment and at the same time, it has made predictions that deviate from GR. Of 
the four possible models, the original version with projectability is not favored by 
observations. On the other hand, the non-projectable version with U{1) symmetry 
proposed by Zhu, Wu, Wang and Shu [85, 84] is free of theoretical inconsistencies 
and compatible with all the observations and in addition, it is embedded into String 
Theory, which was the original purpose of the HL theory. 

Some of the important questions concerning the foundational aspects of the theory 
are worth pursuing. They are the following: 

Renormalizability: Although the theory is shown to be power-counting renormaliz- 
able, establishing full renormalizability has not yet been accomplished. 

Infrared limit: One should rigorously show that HL theory reduces to General Rel¬ 
ativity in the low energy limit using renormalization group flow study. 

Lorentz violating effects: It is interesting to study how the Lorentz violations in the 
gravitational sector propagate to the matter sector and if they have 
any observable effects at low energies. This problem might also help to 
address the question of coupling matter to the theory. 

Here is the list of some important problems that might point us to future directions 
for the development of the theory. 

Quantization: Since HL theory is supposed to be a candidate for quantum gravity, 
it is natural to quantize the theory. This task has already been taken 
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up for some special cases in [147, 148]. It will be very interesting to 
study the graviton-graviton scattering problem and compare it with the 
results obtained in the case of General Relativity by DeWitt. 

Black hole entropy: One obstacle to studying black entropy in HL theory was the lack 
of proper understanding of the black hole horizons, as the theory allows 
superluminal particles that can escape from the Killing horizon. This 
has been remedied with the discovery of the Universal Horizon [149, 
150, 151]. Now the next step is to compute the black hole entropy 
and see how it compares with the standard results of Bekenstein and 
Hawking [152, 153]. 
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